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SUMMARY 

The  goal  of  this  particle  handling  project  was  to  design 
and  construct  an  automatic  particle  sorter.  This  instrument 

was  to  be  capable  of  sorting  out  particles  of  interest  rapid- 
ly and  with  low  loss.  The  particle  diameter  range,  over  which 
this  instrument  was  to  be  effective,  was  5 pm  to  0.1  v*m.  The 
techniques  utilized  were  to  be  highly  independent  of  particle 

snape  or  size  and  present  minimal  damage  and  contamination. 

The  following  techniques  were  investigated  for  use  in  a bread- 
board particle  sorter:  plasma  ashing,  dielectrophoresis,  and 

one-at-a-time  sorting.  In  addition,  particle  aerosolization 
by  plasma  ashing,  gas  dynamic  transport  and  positioning  of 
particles,  and  particle  deposition  from  gas  flow  were  studied 
to  interface  the  three  techniques. 

An  extensive  study  of  the  effects  of  plasma  ashing  on 
particles  and  substrates  was  made.  These  results  should  prove 
useful  in  the  removal  of  particles  from  organic  substrates  and 
in  removing  unwanted  organic  particles  from  samples.  A wide 
range  of  material  and  commercial  substrates  were  ashed  and  a 
large  number  of  gas  compositions  and  relative  ratios  were 
tested.  The  particles  formed  upon  oxygen  ashing  appear  under 
a scanning  electron  microscope  to  be  less  than  0.1  ym  and 
should  be  easily  separated  from  the  larger  inorganic  particles 
remaining.  However,  in  oxygen  plasmas,  samples  of  UO2  showed 
a slight  increase  in  weight  indicating  the  formation  of  higher 
oxides  of  uranixim. 

An  extensive  theoretical  investigation  of  dielectrophoresis 
was  made.  This  resulted  in  the  theoretical  determination  of  par- 
ticle size  dependence,  a dielectric  sorter  module  design,  and 
likely  operating  conditions.  However,  since  the  operating 
conditions  conflicted  with  pressure  and  flow  requirements,  we 
were  not  able  to  test  these  theoretical  results. 

A technique  for  producing  a one-at-a-time  particle  flow  in 
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a gas  was  developed.  Interesting  particles  were  to  be  detected 
by  their  characteristic  fluorescence  which  is  excited  by  a focused 
laser  beam  and  sensed  by  filtered  photo  cells.  The  interesting 
particles  v/ere  then  to  be  deflected  by  an  ion  wind  created  by 
a high  voltage  pulse  between  two  electrodes.  Work  on  the  bread- 
board was  terminated  before  these  ideas  could  be  tested. 

As  advanced  testing  progressed,  it  slowly  became  apparent 
that  the  pressure  and  gas  flow  rate  requirements  of  the  three 
separate  techniques  were  irreconcilably  incompatible.  That  is, 
as  the  likely  operating  conditions  for  each  stage  beccime  known, 
it  became  apparent  that  there  was  no  way  to  interface  them  intb 
one  simple  instrument  as  proposed. 

In  addition,  other  difficulties  came  to  light.  The  dielec- 
tic  constant  to  density  ratios  of  particles  of  interest  were 
found  not  to  be  unique  as  first  proposed.  Plasma  ashing  did  not 
put  the  remaining  particles  evenly  and  compleLely  inco  a flowing 
gas  stream.  As  a consequence  of  these  insoluble  problems,  it 
was  necessary  to  abort  plans  to  construct  a breadboard  instrument. 
For  the  balance  of  the  contract  period,  our  effort  was  directed 
entirely  toward  the  only  stand-alone  component  of  the  system  - 
plasma  ashing. 
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SECTION  I 


INTRODUCTION 


1.0  Hiotoricq],  Perspective 


This  is  the  final  report  for  a two  phase  program  on 
Advanced  Handling  and  Sorting  Techniques  for  Ultrasmall 
Particles".  By  early  1972,  Block  Engineering  had  conceived 
of  a variety  of  particle  handling  techniques  which  seemed 
unique  and  promising.  Those  techniques  were  described  in 
an  unsolicited  proposal'.  That  proposal  led  to  the  funding 
of  the  first  two  phases  of  a proposed  M'ree  phase  program 
to  test  those  techniques,  select  from  among  them  the  most 
useful  ones,  and  construct  a breadboard  automatic  particle 
sorting  instrument. 


In  Phase  I,  a large  number  of  techniques  were  developed 
and  given  preliminary  testing.  The  primary  accomplishments 
are  enumerated: 


(1)  The  diverse  concepts  proposed  were  organized  into 
a unified  system  concept  v;hich  carries  out  the 
desired  functions  in  a logical  order. 


A concept  of  avoiding  liquids  altogether  to 
prevent  contamination  was  developed. 


(3)  Plasma  ashing  was  tested  and  found  to  be  adequate 
for  suspending  the  particles  and  introducing  them 
into  a fluidic  sorter. 


C^)  It  was  demonstrated  that  most  of  the  uninteresting 

.organic  particles  oould  be  eliminated  by  selective 
ashing . 


(5)  A fluidic  system  was  demonstrated  to  be  a low 
loss  technique. 
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(6)  Bulk  sorting  techniques  were  considered  to  eliminate 
uninteresting  particles  early  in  the  fluidic  sorting 
procedure. 

(7)  Suitable  modular  sorter  section  constr-  'ticn  methods 
were  devised. 

(8)  One-at-a-time  sorting  in  a fluidic  system  was  investi- 
gated. 

(9)  Two  nevr  deflection  methods  were  devised  and  tested. 

(10)  Storage  in  the  vapor  state  v/as  demonstrated. 

(11)  Laser  levitation  of  particles  1/10  of  the  diameter 
of  the  smallest  previously-levitated  particles  was 
accomplished. 

(12)  Photophoretic  separation  of  relatively  large 
particles  from  relatively  small  particles  was 
indicated. 

At  the  end  of  Phase  I,  it  was  decided  that  a breadboard 
automatic  particle  sorter  based  on  the  concepts  developed 
during  Phase  I could  be  constructed.  Thus,  the  goal  of  Phase 
II  was  to  construct  a breadboard  particle  sorter  using  the 
following  techniques:  plasma  ashing,  dielectrophoresis,  and 
one-at-a-time  sorting.  In  addition,  the  techniques  of  particle 
aerosolization  by  plasma  ashing,  particle  transport  and  position- 
ing by  flowing  gases,  and  particle  deposition  from  a gas  stream 
would  be  utilized. 

1.1  Performance  Goals  for  Automatic  Particle  Sorter 

A review  of  particle  analysis  problems  determined  for  us 
the  following  specific  requirements  for  this  program  which  have 
guided  our  work. 
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First,  it  is  necessary  that  particle  selection  be  only 
species  dependent  and  be  relatively  independent  of  particle 
shape  or  size. 

Second,  it  is  necessary  that  the  system  have  low  losses. 

To  be  speciJic,  losses  must  be  considerably  fewer  than  losses 
of  available  methods.  It  is  clear  that  currently  losses 
approach  100%  for  particles  below,  say,  0.5  ym.  Thus,  any  sig- 
nificant efficiency  in  the  0.1  ym  to  0.5  ym  range  would  be  help- 
ful. Furthermore,  the  number  of  these  particles  is  likely  to 
greatly  exceed  the  number  in  the  1 ym  to  5 ym  range.  Thus,  al- 
though the  goal  is  no  loss  for  all  sizes,  even  a high  loss  rate 
(e.g.,  20%)  would  offer  a very  dramatic  improvement. 

Third,  damage  and  contamination  must  be  minimized. 

Fourth,  the  sorting  rate  must  be  rapid. 

1.2  Techniques  Chosen  for  Utilization  in  Automatic  Particle 

Sorter 

1.2.1  Plasma  Ashing 

A plasma  is  a state  of  matter  composed  entirely  of  elec- 
trons and  bare  nuclei.  However,  at  very  high  temperatures  or 
in  the  presence  of  strong  electric  or  magnetic  fields,  a glowing 
and  partially  ionized  gas  composed  of  ions,  electrons,  and 
neutral  species  is  also  referred  to  as  a plasma.  In  plasmas 
produced  by  radio  frequency  (RF)  fields,  free  electrons  gain 
energy  from  the  imposed  electric  field  and  lose  this  energy 
through  collisions  with  neutral  gas  molecules.  The  transfer 
of  energy  to  the  molecules  leads  to  the  formation  of  a variety 
of  new  species  including  metastables,  atoms,  free  radicals, 
and  ions.  Thes**  products  are  all  chemically  active  and  can 
be  used  to  react  with  solid  matter  to  form  gaseous  products. 

An  additional,  characterist' c of  such  plasmas  is  the  lack  of 
equilibrium  i etweer  ’ihe  electron  temperature,  T^,  and  the  gas 
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temperature,  T^.  While  the  electrons  are  sufficiently  energetic 
to  cause  the  rupture  of  molecular  bonds,  the  gas  temperature 
may  only  be  slightly  above  ambient  values,  hence  the  term  cold 
or  low  temperature  ashing. 

In  Phase  II,  cold  plasma  ashing  was  to  be  used  in  two 
separate  \/ays.  First,  to  rciiiove  all  the  organic  and  some  of 
the  inorganic  particles  which  contaminate  samples.  While  leav- 
ing particles  of  interest  unaltered,  plasma  ashing  removes  most 
of  the  unwanted  particles  by  reactions  to  form  volatile  products. 
Secondly,  plasma  ashing  appeared  to  be  useful  in  the  removal  of 
the  original  substrate.  As  ashing  progressed,  the  remaining 
particles  could  be  picked  up  by  a flowing  gas  stream  and  car- 
ried into  the  succeeding  stages  of  separation. 

1.2.2  Dio] ectrophoresis 

Dielectric  particles  can  be  deflected  by  inhomogeneous  elec- 
tric fields.  At  low  pressures  (0.1  torr) , this  deflection  is  in 
the  first  order  independent  of  particle  size.  Size  dependence 
enters  as  a second  order  term  v;hich  depends  on  mean  free  path 
and  gas  viscosity.  At  pressui'es  of  10  torr  and  greater,  this 
size  dependent  term  becomes  significant  and  at  higher  pressures 
even  dominant.  In  this  second  stage,  particles  of  one  dielectric 
constant-to-density  ratio  were' to  have  been  separated  from  those 
unwanted  particles  having  very  different  ratios.  After  this 
stage,  those  particles  of  a certain  dielectric  constant  to  density 
ratio  range  will  be  passed  on  to  the  final  stage. 

1.2.3  One-at-a-Time  Sorting 

Once  the  great  bulk  of  uninteresting  particles  have  been  re- 
moved by  the  first  two  stages,  the  number  of  remaining  particles 
should  be  small  enough  to  sort  one  particle  at  a time.  The  par- 
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tides  in  the  transporting  gas  are  so  constricted  that  they 
pass  one-at-a-time  with  constant  velocity  through  a laser 
beam.  The  laser  light  produces  a characteristic  fluorescence 
emission  which  can  be  detected  by  an  array  of  filter  photocells. 
Shortly  downstream, the  interesting  particles  can  be  deflected 
by  a strong  electric  field  into  a separate  gas  stream  for  final 
deposition  or  any  other  type  of  processing. 

1. 3 Organization  of  Phase  II 

The  primary  goal  of  Phase  II  was  to  assemble  and  test  a 
workable  breadboard  particle  sorter  based  on  the  concepts 
examiiicd  in  Phase  I.  The  breadboard  should  have  been  suitable 
as  a b sis  for  a deliverable  model  for  Phase  III.  Phase  II  as 
planr.,cJ  involved  three  separate  steps: 

(1)  To  do  advanced  testing  and  theoretical  investigation 
of  the  three  techniques,  plasma  ashing,  dielectro- 
phoresis, and  one-at-a-time  sorting,  that  are  the 
basis  for  the  three  operating  stages  of  the  auto- 
matic fluidic  sorter.  Also  to  do  advanced  testing 
of  the  interfacing  of  the  three  stages  and  of  the 
final  particle  deposition. 

(2)  To  develop  by  using  the  data  from  step  one,  a de- 
tailed instrument  design  that  is  suitable  for  rapid, 
inexpensive  construction. 

(3)  To  construct  a breadboard  instrument  and  test  its 
capability  and  efficiency  with  standards  and  unknown 
samples. 

Steps  one  and  two  were  to  run  simultaneously  over  the  first 
two  thirds  of  the  contract  period.  These  two  steps  were  to  in- 
teract with  each  other  and  culminate  in  a complete  system  design 
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to  be  built  in  step  three.  The  system  assembly  would  demon- 
strate the  concepts  used  and  the  solutions  to  the  interfacing 
problems.  Of  course,  the  system  would  have  been  improved  con- 
siderably in  transmission  efficiency,  sorting  versatility,  and 
transmission  rate  in  going  from  the  breadboard  model  to  the 
deliverable  system  in  Phase  III.  Nevertheless,  all  of  the  bas- 
ic features  would  have  been  present  in  the  breadboard  model. 

1.4  Results  of  Phase  IT 

An  extensive  study  of  the  effects  of  plasma  ashing  on  par- 
ticles and  substrates  was  made.  The  results  of  this  investiga- 
tion, presented  in  Section  III,  should  prove  useful  in  the  re- 
moval of  particles  from  organic  substrates  and  in  removing  un- 
wanted organic  particles  from  their  samples.  A wide  range  of 
natural  substrates  were  ashed  and  a large  number  of  gas  compo- 
sitions and  relative  ratios  were  tested. 

An  extensive  theoretical  investigation  of  dielectrophor- 
esis, presented  in  Appendix  A,  was  made.  This  led  to  the  deter- 
mination of  particle  size  dependence,  a dielectric  sorter  mo- 
dule design  and  likely  operating  conditions.  Also,  the  problems 
in  utilizing  this  technique  have  been  outlined  in  Section  II. 

As  the  proper  operating  conditions  became  known  for  each 
stage  of  the  sorter , it  became  apparent  that  there  v/as  no  vray  to 
interface  them  into  one  simple  instrirment  as  proposed.  Since 
steps  one  and  two  could  not  yield  a detailed  instrument  design 
suitable  for  construction,  step  three  had  to  be  aborted. 
Consequently,  we  could  not  complete  that  part  of  the  original 
tasking  dealing  with  construction  and  testing  with  both  standards 
and  unknowns  of  a fluidic  particle  sorter.  The  system  design, 
at  the  point  of  termination,  is  presented  in  Section  IV. 
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SECTION  II 

dielectrophoretic  separation 

2.0  Introduction 

The  thorough  theoretical  analysis  of  dirlectrophoresis , 
presented  in  Appendix  A,  is  elucidated  by  a series  of  numerical 
computations  aimed  at  determining  the  theoretical  effectiveness 
of  dielectric  particle  separation  in  our  various  designs.  The 
equations  of  motion  derived  in  Appendix  A are  in  the  form  of 
second  order  differential  equations,  and  are  analytically 
solvable.  However,  th  . approach  taken  here  is  an  iterative 

evolution  method,  well  suited  to  implementation  on  a digital 
computer.  A computer  simulation  of  the  dielectric  separation 
is  presented  with  results  showing  the  separate  effects  of 
particle  size,  initial  position,  and  other  parameters  on  the 

efficiency  of  separation. 

Also,  the  Ijohavior  of  an  ensemble  of  dielectric  particles 
of  two  different  species  (UO2  and  SiO^)  and  of  various  sizes  is 
analyzed.  UO2  and  Si02  were  selected  because  there  was  a great 
difference  in  their  dielectric  constant  to  density  ratios. 
Numerical  results  show  the  theoretically  derived  distribution 
of  these  particles  in  a separation  stage,  and  indicate  the 
degree  to  whic}i  these  species  can  be  separated. 

2 . 1 A Computer  Simulation  of  the  Diclcctrophoreti£ 

Separation  Technique 

2.1.1  Preliminary  Experiments 

In  order  to  gain  insight  into  the  behavior  of  particles 
of  various  sizes,  dielectric  constants  and  densities,  the 
equation  of  motion  for  dielectrophoresis  in  a gas  was  solved. 
The  starting  point  for  these  computations  are  the  equations 
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X = 


--.(K  - 1) 
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(1  - x) 


, 12.3 

{y.  - 2^  ) 
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— 5 ^ T -cd/L. 

pd^  + pdL(a  + be  ) 
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for  the  planar  geometry  and  its  counterpart 


V 

P r,  ^ 


18n  X 


pd^  + pdL(a  -4- 

(2-2) 


for  the  cylindrical  electrode  geometry  (see  Equations  A-39, 
A-40,  and  A-52) . 

Here  x is  defined  as  a normal]  i'.ed  position  coordinate 
(see  Equations  A-37  and  A-38)  and  x and  x are  the  corresponding 
normalized  velocity  and  acceleration. 

Although  Equations  2-1  and  2-2  are  solvable  by  analytic 
means,  it  was  decided  to  simulate  the  behavior  of  a general 
particle  by  computing  its  trajectory  with  a simultaneous 
evaluation  of  the  accelerations  produced  by  the  field  and  by 
the  drag  forces,  in  order  to  decide  v.'lien  drag  forces  became 
considerable. 

The  method  for  tlie  solution  of  the  differential  equations 
2-1  and  2*2  are  very  similar.  In  both  cases  almost  identical 
FORTRAN  source  programs  were  generated.  A listing  of  the 
programs  to  compute  the  trajectories  in  the  planar  and  cylin- 
drical cares,  name.s  FIELD  and  FIELD  3 respectively,  is  attached 
in  Appendix  D. 

In  each  case  a number  of  parameters  are  capable  of  being 
entered  under  program  control,  to  be  used  as  initial  values, 
etc.  These  include  the  initial  position  of  the  particle;  its 
size;  the  applied  voltage;  etc.  Also,  the  geometrical  dimen- 
sions of  the  electrode  configuration  was  variable,  with  the 
inner  electrode  radius  being  fixed  at  1/50  of  the  distance  of 
the  outer  electrode.  Initially,  only  two  choices  for  the 
possible  use  of  the  viscosity  coefficient  were  provided:  one 
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corresponding  to  air  at  standard  temperature  and  pressure,  and 
the  other  corresponding  to  hydrogen  at  1/100  atmospheres  pres- 
sure . 

The  technique  for  tlie  solution  of  the  equations  of 
motion  was  to  update  iteratively  the  position,  velocity,  and 
accelerations  due  to  field  and  drag  respectively,  for  each 
Particle.  For  a convenient  compromise  between  speed  of 
computation  and  the  accuracy  of  the  results,  the  iterations 
wore  done  on  1 nsec  intervals,  v/ith  the  results  being  dis- 
played at  intervale  of  0.1  milliseconds,  or  every  100  iter- 
ations . 

Some  of  the  results  from  executing  the  programs  FIELD  and 
FIELD  3 are  shown  in  Tables  2-1  and  2-5  on  the  following 
pages.  In  each  case,  the  results  have  been  displayed  simul- 
taneously for  identical  particles  under  identical  electrode 
geometry  and  accelerating  voltage  conditions.  In  reading 
these  data,  a minor  discrepancy  botv.'een  the  description  of 
the  equation  of  motion  for  the  planar  geometry  and  its  imple- 
mentation in  FIELD  must  be  mentioned.  In  FIELD,  the  particle 
moves  from  its  initial  normalized  position  (typically  0.95 
cm)  out  to  a norjaalized  final  position  (typically  0.98  cm). 

'ihus  its  velocity  and  acceleration  are  both  positive.  This 
is  in  contrast  with  FIELD  3 where  tlie  particle  is  moving  in- 
wards and  its  velocity  and  acceleration  are  consequently  nega- 
tive. Apart  from  the  minor  differences,  the  absolute  values 
of  the  velocities  and  acceleration  can  be  compared  for  both 
programs  at  each  instant  of  time. 

Tables  2-1,  2-2  and  2-3  show  the  rapidly  increasing 
effect  of  drag  when  the  size  is  reduced,  for  a particle  of 
dielectric  constant  of  20,  and  a density  of  3 gm/cc.  A com- 
parison of  the  field  and  drag  accelerations  indicates  that 
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TABLE  2-1.  Motion  of  a 10  ym  Particle  in  H,  at  1/100  atm 
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TABLE  2-2.  Motion  of  a 5 yni  Particle  in  H~  at  1/100  atm. 
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TABLE  2-3.  Motion  of  a 1 ym  Particle  in  H2  at  1/100 


while  10  pm  particles  are  not  drag-limited  at  1/100  of  an 
atmosphere  H2  pressure,  5 pm  particles  are  more  so,  and  1 pm 
particles  are  almost  completely  governed  by  viscous  drag. 

Table  2-4  indicates  the  effect  of  the  initial  position 
on  the  trajectory  of  identical  particles,  while  Table  2-5 
indicates  its  motion  in  air  under  normal  atmospheric  pressure. 

2.1.2  Determination  of  Particle  Trajectories  for  a Specified 
Distribution  of  Particle  Sizes  and  D ieJectric  Constants 

Following  the  initial  numerical  work,  a more  flexible 
SGfics  of  progranis  was  developed  to  investigate  the  behaviour 
of  dielectric  particle  trajectories.  The  programs  FIELD, 

FIELD  2,  and  FIELD  3,  listed  in  Appendix  C,  each  determine 
different  aspects  of  dielectric  paticle  behaviour.  Each  is 
capable  of  accepting  t)ie  complete  range  of  relevant  variable 
parameters . 

The  new  versioii  of  FIELD  takes  into  account  the  Millikan 
version  of  Stoke 's  Lav:  for  very  small  spheres  in  a gas;  that 
is,  the  mean  free  path  considci'ations  referred  to  in  Section 
A-5.1  have  been  incoj porated  in  the  computation.  Also,  the 
program  is  capable  of  treating  either  the  cylindrical  or  the 
planar  geometry,  with  the  coordinates  being  identically  de- 
fined in  both  cases.  Thus  in  both  modes,  the  particle  is 
travelling  inwards,  towards  x = 0,  and  thus  the  velocity  and 
acceleration  are  both  expressed  as  negative  quantities.  A 
sample  run  to  compare  the  behaviour  of  a lOym  diameter  parti- 
cle in  a)  air  at  standard  pressure  and  b)  in  hydrogen  at 
1/100  atmospheres  pressure  is  shown  in  Table  2-6. 

The  behaviour  of  the  particle  in  air  is  seen  to  be  predomi- 
nantly drag-determined,  whereas  at  the  lower  pressure  in 
hydrogen,  the  drag  effects  are  at  worst  an  order  of  magnitude 
lower  than  field  effects. 
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TABLE  2-4,.  Effect  of  Changed  Initial  Position  of  Particle  Motion. 
1 uin  Particle  in  H_  at  1/100  atm. 
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TABLE  2-6 a Comparison  of  Trajectories  at  Normal  and  Low  Pressures  for  a 10  yim  Particle. 


I The  main  aim  of  this  series  of  simulations  was  imple- 

^ mented  in  the  programs  FIELD  2 and  FIELD  3.  This  aim  consists 

of  generating  the  distribution  of  particles  at  a specific  time 
( instant  as  a function  of  their  size,  dielectric  constant,  and 

density. 

The  program  FIELD  2 is  almost  identical  with  FIELD, 
except  that  the  particle  motion  is  determined  on  a logarith- 
mically varying  time  scale  rather  than  at  equal  intervals  of 
time.  The  minimum  iterative  interval,  as  before,  is  still 
lysec.  Thus  no  additional  error  has  been  introduced.  The 
main  aim  of  this  program  was  to  determine  the  total  time  of 
flight  for  the  particles  suffering  the  highest  acceleration. 
This  then  determined  the  time  for  which  the  entire  range  of 
particles  should  bo  allowed  to  accelerate.  Accordingly,  the 
transit  time  obtained  from  FIELD  2 v/as  entered  into  FIELD  3 
as  a parameter  to  determine  the  distribution  of  tlie  particles . 

' The  program  FIELD  3 determines  the  final  position  for 

particles  of  various  sizes,  densities,  and  dielectric  con- 
stant?;, which  >iave  been  allov;ed  to  accelerate  for  d specific 
time.  At  the  end  of  this  time,  each  particle  position,  velo- 
city, and  accelerations  are  output. 

A typical  output  result  from  FIELD  3 is  shown  on  the 
following  two  pages.  In  these  records,  the  particle  popula- 
tion consists  of  a mixture  of  particle  sizes  ranging  from 
lOym  to  O.lym  in  diameter. 

Also,  two  specific  particle  types  were  selected  for 
i analysis,  namely  particles  of  uranium  dioxide  having  a den- 

I sity  of  10.0  and  a dielectric  constant  of  21.00,  and  particle 

I of  silica,  with  a density  of  2.20  and  dielectric  constant 

1 > 
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2.25.  Note  also  in  these  results  that  in  keeping  with  the 
theoretical  design  criteria  that  were  found  to  be  desirable 
(see  Section  A-7),  the  dimensions  of  the  electrodes  have  been 
diminished  by  a considerable  factor,  with  a total  particle 
transit  distance  being  a maximum  of  100  ym.  The  voltage 
selected  (210V)  is  the  maximum  sustainable  voltage  in  dry 
air  at  one  atmosphere. 

The  result  of  the  outputs  from  FIELD  3 are  summarized 
in  Figure  2-1.  In  this  figure,  the  particles  travel  from 
right  to  loft,  from  the  nominal  distance  of  150  ym,  towards 
the  inner  electrode,  v/hich  is  .located  at  50  ym.  The  shaded 
bar  represents  the  initial  di.stribution  of  particles,  v;hich 
are  depicted  as  starting  with  an  initial  position  uncertainty 
of  -Jt  5 ym.  The  remaining  bars  indicate  the  final  position 
of  the  population  of  particies,  broken  down  by  size,  after  a 
specific  t.lmc  interval.  At  the  end  of  l.)iis  time,  the  popu  • 
lat.Jon  of  10  ym  diaiactcr  UO^  particles  lias  just  reached  tlie 
inner  electrode.  Note  also  that  there  is  an  increasing  spread 
in  the  position  uncertainty  of  each  particle  population  as  it 
travcl.s  tov;ards  the  inner  elci.Lrode. 

Figure  2-1  indicc;te.s  that  , using  the  dimensions  for  the 
electrodes  as  shown  therein,  separation  of  1 ym  diameter 
particles  of  UO2  and  Si02  can  bo  achieved,  provided  the  oper- 
ating pressure  is  1/100  atmospheres  of  II2  and  near  breakdov/n 
voltage  conditions  (in  air)  arc  used.  In  the  presence  of 
submicron  size  particles,  only  incomplete  separations  can  bo 
obtained . 
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J The  species  dependent  factor  in  the  dielectric  acceleration 


equations  (equations  2-1  and  2-2)  is  (K-l)/p.  This  factor  has 
been  calculated  from  the  known  dielectric  constants  and  densi- 
ties for  a nuirJjcr  of  coimnon  inorganic  materials  and  is  tabulated 
in  Table  2-7.  V.’hile  there  is  a large  variance  in  the  dielectric 
constant  (10.7  ± 8.6  or  ± 80%),  there  is  a smaller  variance  in 
our  species  dependent  factor  (2.3  ± 1.3  or  ± 57%).  For  example, 
although  UO2  has  one  of  the  largest  dielectric  constants  (fifth 
largest  on  our  ]ist),  it  stands  33rd  out  of  47  for  the  species 
dependent  factor.  Even  v;orse  is  the  fact  that  it  stands  in  a 
very  dense  portion  of  the  factor  spectrum,  near  many  coinmonly 
occurring  materials.  This  nonuniqueness  of  the  UO2  factor  makes 
its  hu.l  k .sc.parat.i  on  from  other  cormnon  materif'.ls  very  limited  at 
best.  Consequently,  a design  in  v.’hich  there  is  both  a high  and 
a low  cutoff  must  bo  cousiderc;d , ?>ince  UO2  is  neither  near  tlio 
upper  or  lower  limit. 

^ ^ Pre.ssure  Roguirr nvj.nts  for  Effect  ive  ))i electric  Sorting 

Consideriiig  the  cylindrical  geometry,  v;e  again  use  the 
program  FIELD  3 to  detoriiiinc  the  final  position  of  particles 
which  liavo  been  allowed  to  accelerate  for  a specific  time. 

With  helium  a.s  the  carrier  gas  and  a voltage  on  the  central 
wire  of  1C9  volts,  v/hich  is  90%  of  the  breakdown  voltage, 
the  particles  are  brought  in  at  a distance  of  180  ± 5 yra 
from  the  center  electrode.  The  outer  cylinder's  radius  is 
200  yra  and  the  center  wire  has  a 50  ym  radius.  In  Figure 
2-2,  the  final  positions  for  a two  material  system  of  UO2 
(k  = 21,  d =10.9  g/cra^)  and  Si02  (<  = 4.3,  d = 2.6  g/cra^) 
are  shown  for  a range  of  pressures.  At  P = 10  torr,  there 
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TABLE  2-7 


^EPM^TION  RATIOS  FOR  DIELECTRIC  SORTING  OF  COMMOM  matfbtrt c 


I 

1 Material^ 

Dielectric  Constant  (k) 
(at  10®  to  10^  Hz) 

Density  (p) 
(q/cm^) 

(k-1)/P 

I 

1.  Pb(N03)2 

37.7 

4.5 

8.1 

2.  TlCl 

46.9 

7.0 

6.6 

] 

3.  NH^Cl 

7.0 

1.5 

4.0 

n 

1 

4.  *Cassiterite  - SnSnO. 

4 

23.7 

7.0 

3.3 

1 

5 . Na2C03 

8.4 

2.5 

3.0 

6.  *CalcitG  -•  CaCO^ 

8.5 

2.7 

2.8 

■ % 

7.  *Corusnite  - IbCO^ 

18.6 

6.6 

2.7 

« • 

8.  *Tonorite  ••  CuO 

• 

CO 

6.4 

2.7 

* " 

9.  PaCl2 

11.4 

3.9 

2.7 

10.  PbO 

25.9 

9.5 

2.6 

- 

11.  CuSC^ 

10.3 

3.6 

2.6 

A 

12.  NH^Dr 

7.1 

2.4 

2.5 

13.  *Zircon  --  ZrSiO, 

4 

12.0 

4.4 

2.5 

14.  •'■•IlalitG  - Nad 

6.1 

2.2 

2.4 

15.  *7ipalite  - CaF2  ’ 3Ca3p20g 

8.5 

3.2 

2.3 

16.  ’’BaritG  - Ba.SO^ 

11.4 

4.5 

2.3 

17.  PbS 

17.9 

7.5 

2.3 

18.  Tarapacaite  - K2CrO^ 

7.3 

2.7 

2.3 

* 

19.  FeO 

14.2 

5.7 

2.3 

20.  *Dolomite  - CaCO^MgCO^ 

7.5 

2.9 

2.2 

- 

21.  *Anglesite  - PbSO^ 

14.3 

6.2 

2.1 

22.  *Beryl  - 3Be0’Al203 • 6SiO 

2 6*6 

2.8 

2.0 

23.  *Gypsum  - CaS0^*2H20 

5.7 

2.3 

2.0 
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TABLE  2-7  (continued) 


I • 


1 24.  AI2O3  - a 

0.8 

4.0 

2.0 

|25.  CaP2 

7.4 

3.2 

2.0 

26.  *Sylvite  - KCl 

5.0 

2.0 

2.0 

27.  *SaltpGter  - KNO, 
1 J 

5.0 

2.1 

1.9 

,28.  KjCOj 

5.6 

2.4 

1.9 

j 29.  *Smii:h5onitG  - ZnCO^ 

9.3 

4.4 

1.9 

i30.  KClO. 
1 

5.1 

2.3 

1.8 

31.  *Arconite  - K,,SO. 

i 

5.9 

2.7 

. 1.8 

|32.  *CGrargyrite  ••  AgCl 

11.2 

5.6 

i.si 

'33.  UO2 

I 

21.0 

10.9 

1.8 

34 . >’-Aj.'acjonitG  - CaCO^ 

6.1 

2.9 

1.8  1 

35. 

5.2 

2.3 

1.8  1 

j 36.  .i.i.G  • 7nj}'.r 

12.2 

• 

6.5 

1.7 

37.  *Ka  Unite  - l'J\l  (SO^ ) 2 • I2H2O 

3.8 

1.8 

1.6' 

|38.  'Ha.l aca’iitc  ••  CnCO^Cu  (Oil)  2 

7.2 

4.0 

1.6  ^ 

39.  ’'Alrrn.nin  - Al„0-, 

6.5 

3.7 

1.5 1 

40.  Ha{N03)2 

5.9 

3.2 

1.5  i 

j41.  *Quartz  - Si02 

4.3 

2.6 

1 

1 

1.3 

'42.  *Diojnond  - C 

5.5 

3.3 

1.3  * 

.43.  KI 

1 

5.6 

3.6 

1.3 

|44.  AgCN 

5.6 

4.0 

j 

1.2 

45.  Hg2Cl2 

9.4 

1.2 

1.2 

46.  *Sand  - Si02 

2.6 

2.6 

0.62 

47.  *Cotunnitq  - PbCl2 

4.2 

5.9 

0.55| 

+ Chemical  names  indicate  data  for  very  pure  raaterials 
Minerals  and  other  naturally  formed  materials 
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is  a partial  separation  like  that  shown  in  Section  2.1. 

However,  this  separation  is  too  little  for  any  practical  use. 

At  P = 1 torr,  there  is  sufficient  separation  for  a two 
material  system.  However,  only  at  P = 0.1  torr  is  there 
sufficient  separation  for  a multimaterial  system. 

In  Figure  2-3,  we  show  the  results  of  a multiparticle  separa- 
tion for  1 ym  diameter  particles  with  an  applied  pressure  of 
0.1  torr.  If  the  particles  falling  betv;een  90  ym  and  120  ym 
from  the  center  of  the  wire  were  -eparated  from  the  remainder 
of  the  particles,  there  is  useful  separation  of  UO^  from  par- 
ticles having  a (k-1)/p  factor  greater  than  2.3  or  less  than  1.4 


In  conclusion,  these  results  show  that  an  effective 
separation  can  only  be  possible  at  pressures  less  than  about 
0.1  torr  (100  micron  of  Hg) . 

2 . 5 lKabora^o2:v_J’>^.Pi'.ri^^^ 


A simple  qualitative  c>:periit*ent  v;as  conceived  to  demon- 
strate dielectric  deflection  of  a slowly  moving  beam  of 
particles.  A very  fine  particle  stream  of  lov;,  but  unmeasured, 
particle  density  was  generated.  The  followii.j  equations  from 
Appendi::  A are  used  to  calculate  the  parameters  of  the  e::peri- 

ment. 


Co(K-l) 


o_ 

P'r. 


(2-3) 


where 


P'  = £n(rj^/r^) 

with  r = 0.24  cm  (inner  electrode  radius)  and  r^^  = 2.54  cm 
(outer  electrode  radius) , an  applied  voltage  (V^)  greater 
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than  5.7  kV  should  deflect  UOj  particles  moving  at  a velocity, 
Vp,  of  about  3 cm/sec  tangential  to  the  inner  electrode  (see 
Figure  2-4).  Operating  the  system  at  7 kV  and  Vp  = 3 cm/sec, 
extensive  deflection  of  a beam  of  1 pm  UO2  particles  was  ob- 
served. It  was  not  possible  to  test  the  particles  for  surface 
changes.  On  a filter  paper  used  to  collect  the  particles,  it 
was  observed  that  UO2  particles  deflected  rather  evenly  with 
a tapering  off  starting  2 cm  from  the  point  of  no  deflection. 
The  same  experiment  was  then  tried  using  1 pm  Si02  particles. 
With  Vp  again  at  = 3 cm/sec,  there  was  little  deflection  of 
the  Si02  until  voltage  was  increased  to  10  kV.  The  quali- 
tative results  of  these  tests  are  shown  in  Table  2-8. 

TABLE  2-8 

QUALITATIVE  ItBSULTS  OF  DIELECTRIC  EXPEBIM-ENT 


Voltage  (kV) 

1 ym  UO2 

Er.1.,1  mote  of  the  A-mouni 

1 ym  Si02 
of.  Def leci-.ion 

4 

None 

None 

5 

None 

None 

6 

Some 

- 

7 

Some 

None 

8 

Extensive 

None 

9 

Extensive 

Some 

10 

- 

Some 

11 

— 

Extensive 

Since  at  8 kV,  there  appeared  to  be  extensive  deflection  of 

UO2  and  almost  none  for  Si02,  a separation  experiment  was 

performed.  A mixture  of  UO2  and  Si02  was  placed  in  the  particle 

chamber.  With  V = 8 kV,  a 1:1  mixture  of  1 ym  UO,  and  1 ym 
o ‘• 
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SiO^  particles  showed  slight  deflection  and  no  separation. 

Since  the  particles  under  a microscope  appeared  to  be  highly 
conglomerated,  a second  particle  generator  was  constructed 
so  that  the  UO2  and  EiO^  could  be  combined  in  the  low  density 
gas/particle  beams.  The  results  were  nearly  identical  to 
those  usijig  the  intimate  mixture. 

Under  a very  intense  lamp  the  individual  particles  could 
be  observed.  These  particles  demonstrated  no  acceleration 
toward  the  center  high  voltage  electrode.  Even  at  15  kv, 
no  migration,  either  dielectric  or  Coulombio,  was  ob.served  by 
tho  particles  moving  very  near  the  wire.  Although  microscopic 
examination  of  the  particle.s  showed  that  they  were  still, 
although  less,  conglomerated,  their  central  acceleration ’should 
not  have  been  effected.  Consequently,  the  viscous  drag  in 

nitrogen  and  in  helium  (the  two  gases  used  in  these  experiments) 

r.t  atnoepherne  pror:."ur''  i-.  •f-.-n  

f ujTs,  i.u  too  xor  us  -Lo  observe  any 

<3  i p 1 c c 1-  r i c (j  o f .1  c c I-.  i r. : 1 , 


The  deflection  that  we  observed  in  our  first  experiments 
wore  due  to  t)ie  gas  dynamics  in  our  enc.losed  chamber.  The  dif- 
ferent results  with  our  two  mater.lals  were  probably  a result  of 
their  varying  densities. 


From  these  experiments  v?o  can  conclude  two  things: 

(a)  dielectric  deflection  can,  if  ever,  be  observed 
only  at  very  low  pressures,  and 

(b)  our  calculation  of  dielectric  forces  may  be  greater 
than  they  actually  are. 


Dielectric  Forces  versus  Coulombic  Force.g 

In  Appendix  A,  we  showed  that  the  dielectric  force  is 

^d  = E(3E/8r)V  (2-5) 

where  V is  the  particle  volume. 
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It  is  well  known  that  the  Coulombic  force  is 


Fc  = Eq 


(2-6) 


For  the  Coulombic  force  to  be  of  the  same  order  magnitude  as 
the  dielectric  force  (i.e.,  F^/f^  = 1)  then 

q = no  = Cq(k-I) (3E/3r)V 

where  n is  the  number  of  electrons  whose  charge  is  e.  Assuming 
spherical  particles  and  using  the  cylindrical  geometry  and 
equation  A-33, 


r An(r3_/r^) 


(2-7) 


with  V - 7Td  /6,  the  number  of  eloct.vons  that  will  give  a 
Coulombic  force  equal  to  the  dielectric  force  is 


,3 

2 " ■ ^ 
6er  iln(rj|/r^) 


(2-8) 


bet  V - 190  volts,  r - 175  yra,  r^  50  ym  and  r^^  - 200  ym, 
v;hci  for  UO2,  < ^ 21.  Table  2-9  shows  the  number  of  electrons 

TABLE  2-9 

THE_  NUMBER  JjigCTOpjjS^Tn^^^  A COiyLO>mj^FQRCE 
EQUAL  TO  THE  DIELFCTRIC  FORCE 


Particle  Diameter  (um) 

50 

10 

5 

3 

1 

0.8 

0.5 

0.1 


n (^  of  electrons) 

320,000 

2,600 

320 

70 

2.6 

1.3 

0.32 

0.0026 
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for  particles  of  UO2  of  various  sizes.  For  particles  of 
diameter  less  than  1 um,  any  surface  charge  will  give  a large 
Coulombic  force. 

Particles  in  the  atmosphere  are  known  to  be  charged; 
electrostatic  precipitators  work  very  effectively  on  the 
principle  that  all  particles  have  some  charge.  Charge  is 
caused  by : 

(a)  Friction  - particles  colliding  with  each  other  and 
with  gas  molecules. 

(b)  Reactions  - cliemical  and  nuclear 

(c)  External  forces  - ultraviolet  radiation,  gamma 
radiation,  ionic  discharge. 

Even  if  the  particles,  which  v;e  must  assiime  are  likely  to 
be  charged  v;hen  they  ieave  the  plasma  asher  section,  could  be 
ef f C-Cti.vei y discharged,  v/ithout  losing  them  on  the  discharging 
surfaces,  it  is  extremely  vinlikely  tJ.ey  could  be  neutralized 
and  kejjt  neutral  to  the  degree  v.’here  there  v/ould  literally  be 
no  free  charges  on  their  surface. 

A solution  v.’cuid  be  to  rapidly  oscillate  the  electric 
field.  In  order  to  see  a relatively  uniforra  field,  the  field 
roust  oscillate  at  least  once  each  time  a particle  trans- 
verses  its  diameter.  For  1 pm  particles  moving  at  1800  cm/sec, 
the  freguency  must  be  1.8  x 10^  Hz.  This  rapidly  oscillating 
field  creates  a situation  similar  to  our  RF  plasma  stage.  If 
the  pressure  is  low,  i.e.,  below  10  torr,  a plasma  discharge 
will  occur.  Consequently,  this  consideration  dictates  a 
higher  pressure. 
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SECTION  III 
PLASMA  ASHING 


3.0  Introduction 


In  an  RF  oxygen  plasma,  the  primary  reactive  species  has 
been  determined^  to  be  0^  molecules  and  0 ^P  oxygen  atoms, 
their  individual  concentrations  ranging  from  10  to  20%  depend 
on  the  dischaige  conditions.  Both  the  atomic  and  excited 
mclecular  oxy(}en  react  with  organic  substrates  and  gas  phase 
species.  In  the  glow  region  of  the  plasma,  there  is  a free 
electron  concentration  of  approximately  10^^  cn“^.  However,  a 
few  millimeters  away  from  this  glow  region,  nearly  all  of  the 
ionic  species  have  disappeared  but  the  other  "active"  species 
persist  for  a considerable  distance  before  they  are  deactivated 
either  reactively  or  by  collision. 


The  direct  generation  of  oxygen  atoms  in  the  discharge 
glo'.i’  region  rcr.u3tr.  in  part  from  o3.cctrcn~molcculc  reactions, 
molecule  ion- electron  reactions,  and  electronically  excited 
molecular  oxygen  dissocation.  High  concentrations  of  atomic 
oxygen  can  be  attained  by  the  catalytic  effects  of  other 
foj eign"  gas,  such  as  hydrogen,  nitrogen,  or  water  vapor. 
Those  impurities  provide  many  of  the  kinetic  pathways  leading 
to  atomic  oxygen  production.  Nitrogen  as  an  impurity  in 
discharged  oxygen  gives  such  reactions  as: 


NO^  + e~  N + 0 
N + O2'’’  no'*’  + 0 


(3-1) 


In  a very  pure  oxygen  discharge,  there  is  the  usual  glow,  but 
the  plasma  v/ill  contain  a very  low  level  percentage  of  oxygen 
atoms . 


i 


1.  J.  R.  Hollahan  and  A.  T.  Bell,  Techniques  and  Applications 
of  Plasma  Chemistry,  John  Wiley  & Sons,  Inc.,  1974. 
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3 . 1 Optimum  Conditions  f o r P.l at^ma  Ashing 

3.1.1  Pressure  and  Flov;  Rato 


A fluidic  system  v;as  developed  v.’hore  tlio  pressure,  flow 
rate,  and  gas  mixture  could  bo  independently  controlled.  This 
system  is  shown  in  Figure  3-1.  The  aslicr  is  a Togal  PLASMOD 
whoso  RF  pov.'or  is  continuously  variable  from  0-100  at  13.56 
MHz.  The  PLASKOD  was  fitted  with  3.25"  ID  by  6.25"  long  Pyrex 
reaction  chamber.  The  meclianical  pump  is  a Welch  14  02H  wlii  ch 
has  a pumping  speed  of  160  liters  per  minute.  The  pressure 
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TIjc 

fol 

low 

ing 

c , s : 

f.r'rtur.', 

we;’- 

e .'  tudied 

1 

(a) 

100?. 

^2 

(b) 

96? 

°2' 

4%  1 

(c) 

00? 

20% 

(d) 

50% 

50% 

N2 

(o) 

48% 

' 

40% 

N2 

, 4%  1 

CF, 

(f) 

50% 

®2' 

50% 

»2 

For  each  mixture,  two  samples  each  of  approximately  2 grams 
of  Whatman  #44  (Ashless)  filter  paper  were  used.  The  results 
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PLASMA  ASHER  AND  FLOW  SYSTEM 
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are  chown  in  Table  3-1  and  in  Figure  3-2.  In  all  canec,  the 
variation  in  the  percent  residue  after  each  time  interval  was 
within  ± 0.05%  of  the  average.  This  was  small  compared  to  the 
relative  differences  between  mixtures.  The  1:1  mixture  of  O 
and  ^2  v;ith  CF^  v;as  the  fastest  mixture  tested.  The  addition 
of  4%  CF^  did  increase  the  rate  slightly  but  insignificantly. 
Adc  1 ng  nitrogen  to  oxygen  markedly  increased  the  rate  of  ash- 


TAEJ.F  3-j^ 

AS)1 1 NG ItM’]'  I'0)<  V/'RTQUS  G7  S ’ ' IX'J'UFJJS 


1 

Ca5 

Mi  xl.urc 

^2 

^2 

\7  CJ\ 

'i 

4 : 1 

1:1 

1:1 

1:1 

Average 

2 Hour  1 

58. 

32 

52. 

94 

38  . 

75 

34  . 

81 

32. 

17 

42. 

0 3 

4 Hour  6 Hour 
33.68  20.51 


2 8.4  r 
15.1.0 
11.87 
9.8C 
18.^2 


3 5.03 
5.61 
3.  81 
3.1.1 
7.92 


8 Hour 

10 

10.18 

4. 

7.32 

3. 

1.4  2 

0.. 

1.35 

0.: 

kT 

CO 

o 

o.c 

1 2.34 

1 .] 

A more  important  tost  is  the  relative  completeness  of 
ashing,  as  well  as  the  inertness  of  UO.  to  the  mixture.  Sample 
of  cotton  cloth  (CC)  and  UO^  were  subjected  to  plasmas  of  the 
follov.’ing  mixtures  for  18  to  24  hours: 


(a) 

100%  O2 

(b) 

96% 

®2' 

4% 

CF4 

(c) 

50% 

^2' 

50% 

N2 

W) 

CO 

O2, 

48% 

N2 

(e) 

effi 

0 

CO 

®2' 

20% 

CF 
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O results  arc  shown  is  Table  3-2.  The  addition  of  CF,  gives 
consrderably  loss  ash  fro.n  eotton  cloth.  While  the  l = /o'/„ 
grves  a ™ch  taster  ashing  than  pure  0,,  the  final  percen  ^ 
-srdue  rs  nearly  the  sa.e.  The  1:1  with  CF,  gives 

^th  the  fastest  and  most  complete  ashing,  however,  tte  test 
for  corrosron  to  ho,  shows  that  mixtures  containing  CF,  caus^ 
a rerrous  loss  in  weight  to  UO,  samples,  since  plasma^  can 
ly  etch  away  at  the  surface,  a .small  weight  loss  in  this 
experrment  indicates  that  an  exposed  particle  10  urn  or  less 
rn  drameter  would  be  completely  destroyed  if  CF,  in  any  con- 
ten  ration  as  present  in  the  plasma  gas.  Consequently,  CF 
should  not  be  used  on  samples  where  low  losses  are  reguirel 

Samples  exposed  only  to  pure  0^  or  to  O2/N.  mixtures 
a owed  a slight  increased  weight.  This  is  undoubtedly  due  to 
e formatron  of  higher  o.xides  lihe  UQj.  since  no  materia]  was 
lost,  the  1:1  o„/K2  mixtures  wore  deemed  the  best  to  use  in 
the  rcnininder  of  our  oxper i rn^'fital  v;ork. 


3.1.3  n>;perimrnt  al  Procedure 

The  samples  were  a.shcd  in  Pyrox  petri  dishes  60  mm  OD  x 
15  nmi  height.  Before  initial  weighing,  the  dishes  were  placed 
an  the  operating  asher  for  at  least  1 hour  to  remove  all  ad- 
sorbed water.  AIJ  weighing  was  do,)c  on  a Mettler  H15  analyti- 
cal balance  which  has  a sensitivity  of  0.1  mg  and  a precision 
(standard  deviation)  + 0.05  mg. 

All  samples  were  dried  for  at  least  24  hours  in  a des- 
iccating cabinet  in  which  the  desiccant  drierite  was  used. 

A silica  gel  desiccant  was  placed  in  the  balance  chamber  to 
reduce  moisture  adsorbtion  during  weighing. 
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TABLE  3-2 

PERCENT  RESIDUE  AFTER  LONG  TEIM  ASHING 


Gas 

°2 

V CF, 

1:1 

1:1 

02/11, 

1 : 1 

v; 

C^4 

1:1 

Oj/Hj 

1:1 

0^/;j2 

cr, 

°2 

°2  CF4 

1:1 

1:1 

"z/'T 

1:1 

O2A2 

W 

CT, 

1:1 

O2/N2 

w 

CF4 

1:1 

O2/N2 

w 

CF, 

4;  1 

Oj/CF^ 

4:1 

Oj/CF^ 

Ashing  Time  Initial  v;t. 
Sample*  hours  gm 


23. G 
22.1 
21.4 
21.4 

24.0 

18.1 
iC.l 
23.  C 
22.1 
21.3 
21.3 
24 . 0 
20.2 
20.2 
21.3 
21.3 


3. 58G3 
2.8164 
2.1107 
l.lSf)9 
0.7992 
3. 9296 
3.9696 
8./' 692 
7.9136 
6.6439 
7.3126 
9.1268 
13.1687 
7.8642 
6.3874 
10.1896 


Fj.nal  V7t 


0.0114 
0.0042 
0.0061 
0.0028 
0.0005 
0.0016 
0. 00.'’8 
8. 5810 
7.3090 
6.7250 
7.3777 


11.4950 
7.1847 
4 . 6238 
7.6249 


% Residue 


0.32 

0.15 

0.29 

0.24 

0.06 

0.04 

0.07 

101.32 

92.36 
101.22 
100.89 

89.36 
87.29 

91.36 
72.39 
74.83 
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Th?  UO2  samples  v;tre  heated  in  an  oven  at  150“C  for  at 
least  24  hours  before  ashing.  The  UO2,  99.8%  pure,  was 

purchased  from  Research  Organic/Inorganic  Chemical  Corporation, 
Belleville,  NJ  07109. 

During  each  run,  only  two  samples  and  a blank  di jh  could 

be  placed  in  the  asher.  Table  3-3  shows  the  initial  and  final 

weights  of  blanks  for  several  runs  using  different  gases.  For 

those  runs  in  which  CF^  was  present,  there  was  an  average  weight 

loss  of  -1.9  mg.  The  sample  weights  in  Table  3-2  v;ere  corrected 

for  a -1.9  mg  loss  in  v/eiglit  by  the  container  in  those  runs 

where  CF^  was  present.  In  the  other  runs,  there  was  a standard 

deviation  ± O.IG  mg.  The  maximum  error  in  estiinated  residue  i- 
± 0.016%. 


TABLE  3-3 


VARI7\TIO?;  T TG3IT  OF  BLANK S 


Ga  fj 

Ashing  Time 
hours 

1 liiitral  WtT* 
gm 

Fina'l  Wt. 
gm 

A (mg ) 

°2 

23.6 

16.1312 

16.1313 

+0.1 

O2  w CF4 

22.1 

17.1092 

17.1074 

-1.8 

1:1  O2/N2 

12.8 

16.1312  • 

16.1312 

0.0 

1:1  O2/N2 

21.4 

16.1314 

16.1314 

0.0 

1:1  O2/N2 

18.8 

11.5118 

11.5115 

-0.3 

1:1  O2/N2 

20.3 

16.7048 

16.7049 

+0.1 

1:1  O2/N2 

19.8 

15.6736 

15.6738 

+ 0.2 

1:1 

17.0 

14.4265 

14.4267 

+0.2 

1;-1  O2/N2 

17.7 

17.5953 

17.5953 

0.0 

1:1  O2/N2  w CF^ 

24.0 

16.1315 

16.1298 

-1.7 

1:1  O2/N2  w CF^ 

18.1 

16.7894 

16.7875 

-1.9 

1:1  O2/N2  w CF^ 

19.2 

16.9765 

16.9742 

-2.3 
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In  all  runs,  the  plasma  was  tuned  to  maximum  power  ab- 
sorbed. The  characteristic  plasma  colors  of  the  following  pure 
gases  are; 

Oxygen  Pink 

Nitrogen  Orange 

Carbon  Dioxide  Blue 

Hydrogen  Redish  Pink 

Carbon  Tctraf luori do  Bluish  White 

V7hcn  an  organic  sample  begins  ashing  in  an  oxygen  or  oxygen/ 
nitrogen  plasma,  the  color  is  initi^illy  blue  due  to  the  CO2 
formed.  However,  this  .shortly  (10  to  20  minutes)  changes  to 
tlie  char ac tor i.s tic  pink  of  oxygen.  Consequently,  color  has 
no  utility  in  determining  the  complc?teness  of  ashing.  It  is 
useful  in  dotcriai ning  .if  the  asher  is  well  tuned. 

3.2  Oxyf’on  of  n ‘ e_r;"i  eri  .-O  s 

3.2.1  Preliminary  Experiments 

Tabac  3-4  .sho\7s  tlio  results  of  sJiort  t-shiJigs  (8  to  10  hrs.) 
v;ith  pure  O^.  All  runs  made  under  tliese  eonditions  arc  listed. 

The  prc.ssurc  was  approxiinat.c.ly  3 torr  and  the  surface  temperature 
was  estiraated  to  be  about  120°  C.  These  experiments,  not  done 
with  tlie  systcra  as  sliown  in  Figure  3-1,  had  the  gas  going  directly 
through  a needle  va.ivc  to  the  reaction  chainJaer.  Consequently, 
the  flo.;  rate  was  not  monitored. 

3.2.2  Later  Experiments 

In  the  experiments  shown  in  Table  3-5,  the  system  was  oper- 
ated as  described  in  Section  3.1.1  with  a 1:1  O2/N2  mixture. 

For  those  ashings  which  were  visibly  incomplete,  the  final  weight 
is  marked  by  an  I.  All  runs  made  have  been  listed  and  there  i.s 
a very  good  agreement  between  the  results  of  identical  runs. 

When  two  and  three  samples  were  run,  the  observable  rate  of  oxi- 
dation varied  indicating  that  the  power  was  not  ex'enly  distri- 
buted throughout  the  cavity. 
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TABLE  3-4 


ASHING  OF  ORGANIC  I4ATFRIALS  V7ITH  OXYGEN 


o T w . . T Initial  Wt. 

Sample  Material 

^ gm 

Final  Wt. 
gm 

% Residue 

A.  Pure  Sub.stanccs 
1.  Polyotliylenc 

1.0022 

0. 0012 

0.1 

2. 

Poly viny Ich ] oride 

0.5029 

0.0031 

0.6 

3. 

Ny  Ion 

1.5160 

0.0025 

0.2 

4. 

Graphite 

1.0015 

0.  0000 

0.0 

5. 

S tarcli 

1.0244 

0.0002 

0.0 

B . Comp  lex  Ma  to  r i a .1  s 
1.  Huiiu'ui  Hair 

1. OOCO 

C.045;: 

4.9 

2. 

Pine  Needier; 

1.983G 

0.0719 

3.G 

3. 

Alldlfa  (dry  pov.'der) 

2.0006 

0.1445 

7.2 

4. 

Maple  Leaf 

0. 5901 

0.0724 

12.3 

5. 

Boo);  Paper 

G. 3129 

0.0892 

1.4 

G. 

Imping  Paper 

1.031G 

0.0089 

0 . 3 

7. 

Filter  Paper  (Cotton) 

2.7250 

0.0089 

0.3 

8. 

Coal  Dust 

0 . 5005 

0.0392 

7.8 

9. 

Pov.’dc.rf  ('  VJhole  Milk 

1.0013 

0.0657 

6.6 

10 

Motor  Oil  (Clean) 

1.0003 

0.0003 

0.0 

11. 

Sav;dust 

1.0048 

0.0304 

3.0 

TABLE  3-5 


1 

ASHING  OF  ORGANIC  MATERIALS  WITH  3:1 

OXYGEN/NITROGEN  MIXTURE 

1 

f 

Sainjilo 

Material 

Asliincj  Tii.ie 
hours 

Initial  Wt. 
gm 

Final  Wt. 
gm 

% Residue 

! 

1 

A.  Filter  Papers 

1. 

Whatman  ^50 
(Hardened--0 . 025?,) 

49.5 

13.2553 

0.0089 

0.067 

1 

2. 

Whatman  i;541 
(Hardened  Ashless- 
O.OOS?,) 

49.5 

8. 6G90 

0.0002 

0.002 

1 

3. 

Whatiaan  tM4 
(iAr,liless-0 . 03  v) 

49.5 

9.5029 

0.  OOIG 

0.017 

12.8 

4.3740 

0.0017 (I) 

0.039 

f 

Mil  3 i pore  Duralon 
(NC-Nylon) 

18.8 

0.2008 

0.0002 

0.10 

5. 

Millipore  Ce3 elate 

•.r 

i 

G. 

(E7i"Ccl l.ulose  /acetate) 
Mi]  l.i.porc'  Polyvic 

18.0 

0. 20G2 

0.0004 

0.19 

- 

(BS-l’e/lyvi  ny3. 
Chloride ; 

1 8 . 8 

0.3338 

0.0664 

TO  n 

7. 

Mr-  Mid  l.-ipore 
( S !i ■ i -1  i : e d E s c r of 
Cellulose) 

12.  C 

0.2240 

0.0001 

o 

o 

13.  Clcj Hi jjifj  ’•ii<  i'.or  i a.l 

1. 

lOO'o  Cotton 

21.4 

2.3107 

0.0061 

0.29 

21.4 

1.3  5 59 

0.0028 

0.24 

2. 

75?.  ]<:,yo;i/25t  Cotton 

20.3 

3.111.4 

0. 0167 

0.54 

17.7 

2.465G 

0.0134 

0.54 

3. 

50?,  Cotton/50?, 
Polyester 

19.8 

1.323  2 

0.0246(1) 

1.86 

17.0 

0.8214 

0.0531  (I) 

6.46 

1 

17.7 

0.9836 

0.1414  (I) 

14.37 

1 

4. 

50i  Orion  Acrylic 
25%  Stretch  Nylon 
25%  Cotton 

13.8 

1.9510 

1.0083 

0.43 

17.0 

1.5008 

0.0062 

0.41 

13.8 

2.5006 

0.1511(1) 

6.16 

5. 

100%  Orion 

17.0 

1.1G97 

0.0104 

0.89 

20.3 

2.2564 

0.0201 

0.89 

6. 

Irish  Linen  (flax) 

17.7 

1.6122 

0.0042 

0.26 

(* 

19.8 

1.3185 

0.0035 

0.27 

i 


iii-ii 


Block  l.n/tinrrrinjt,  Inc. 


TABLE  3-5  (continued) 


7. 

60%  V/ool,  5%  Cotton 
5%  Rayon,  10%  Acetate 
20%  Stretch  Nylon 

35.6 

35.6 

8. 

Icelandic  Wool 

21.4 

23.1 

9. 

Australian  Wool 

21.4 

23.1 

Natural  Matorials 

1. 

Bark 

24.5 

48.3 

2. 

Dried  Jjoaves 

7.1 

24.5 

3. 

Frcsli  Loii.vec 

22.1 

22.1 

Tins 

ues 

1. 

Kirav.’ipen 

20. G 

48.3 

2. 

Klc'onr':: 

20.6 
4 8.3 

2.6771 

2.1087 

0.0461 

0.0224 

1.72 

1.06 

3.0476 

1.6308 

0.0290 

0.0145 

0.95 

0.89 

2.7511 

2.2460 

0.0167 

0.0124 

0.61 

0.55 

0.8829 

2.3694 

0.0529  (I) 
0.0242 

5.99 

1.02 

0.0855 

0.6831 

0.0043 

0.0092 

5.03 

1.35 

0.7836 

0.8643 

0.0061 

0.0084 

0.78 

0.97 

1.0871 

1.1242 

0.0257  (I) 
0.0121 

2.36 

1.08 

0.9718 

1.4859 

0.008- 

0.0131 

0.91 

0.88 
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3.2.3  Photographic  Exo.mination  of  the  Residues 

Although  the  residues  v;ere  not  chemically  analyzed,  they 
are  typically  oxides  of  trace  elements  such  as  Ca , Cu,  Fe,  Mg, 

Zn,  Ma,  and  K^.  The  particles  formed  upon  ashing  are  typically 
considerably  smaller  than  0.1  pm.  Figure  3-2  shows  430X  en- 
largement of  ashed  and  crushed  sample  Millipore  Duralon.  Since 
no  large  particulate  matter  can  be  identified  and  since  a 1 ym 
particle  would  give  a .43  mm  image  on  our  film,  this  picture  shov;s 
that  particles  formed  are  well  less  than  1 Figure  3-3  shows 

a second  sample  of  ashed  Duralon  which  liad  10  ym  and  smaller 
UO2  particles  on  its  surface.  These  particles  show  up  clearly 
under  a microscope.  Figures  3-4  and  3*5  show  that  the  plasma 
does  not  change  the  material  matrix  in  Irish  linen,  and  Figures 
3-6  and  3-7  show  the  same  effect  for  cotton  cloth.  Figures  3-8 
and  3-8  show  lOOt  orlon  acrylr'c  before  and  after  ashing.  Fig- 
urers  3-10  and  3-13  shov;  Australian  v’ool  before  and  after  ashing. 


3.3 


the 

all 


Ashing  of Mater .ials 

Table  3-G  shov;s  the  result  of  preliminary  as}xings  under 
same  conditions  as  described  in  Section  3.2.1,  Data  from 
runs  made  tire  shov/n. 


1.  J.R.  Hollahan  and  A.T.  Bell,  Techniques  and  Applications  of 
Plasma  Chemistry,  John  Wiley  & Sons,  Inc.,  1974. 
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FIGURE  3-3 

[ ASHED  ANO  CRUSHED  DURLON  4 3 OX 


FIGURE  3-4 

ASHED  AND  CRUSHED  DURLON  WITH  10  pm  PARTICLES  4 3 OX 


I 


I 
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FIGURE  3-6 

ASHED  IRISH  LINEN  lOOX 
IIi-15 
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FIGURE  3-7 
COTTON  CLOTH  lOOX 


FIGURE  3-8 

ASHED  COTTON  CLOTH  lOOX 
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FIGURE  3-10 

100%  ORLON  ACRYLIC  lOOX 
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figure  3-11 
AUSTRALIAN  VJOQL  10 OX 


i ( 

I 

I ' FIGURE  3-12 

I ~ ' 

( ASHED  AUSTRALIAN  WOOL  lOQX 

II  i-18 
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TABLE  3-6 


ASHING  OF  INORGANIC  MATERIALS  WITH 
OXYGEN,  HYDROGEN  AND  TETRAFLUOROMETHANE 


Sample 

Material 

Initial  Wt. 
gm 

Final  Wt. 
gm 

% Residue 

A.  Ashed  v;ith  0>rygen 

• 

1. 

Clay  Encrusted  Sand  1.2715 

1.1845 

93.2 

2. 

Potting  Soil 

5. 5802 

4.0943 

73.4 

3. 

Silica  Dust  (SiO^) 

0. 5609 

0.5601 

99.9 

4. 

Zinc  Oxide 

4.1525 

4.1506 

100.0 

5. 

Zirconium  Nitride 

0. 5609 

0.3409 

60.8 

6 • 

Uranium  Dioxide 

2.9836 

2.9834 

100.0 

B.  Ashed  with  Hydrogen 

1. 

Silica  Dust 

0.6304 

0. 6 301* 

100.0 

2. 

Uranium  Dioxide 

3.1567 

3.1564* 

100.0 

Ashing  times  v?ere  about 

18  hours. 

C . As 

hed  v.'ith  Tetraf luoromethane  (CF^) 

1. 

Sii.ica  Dust 

0.5173 

o.oooo"^ 

0.0 

2. 

Uranium  Dioxide 

1. 4142 

1.0960’^ 

77.5 

I j 


Adjusted  for  loss  in  v;eight  of  the  pyrex  container.  Ashing 
times  here  v;ere  1 hour-. 


1 

■I 

'I 


I 

I 

I 3.4  Cone] us ions 

Plasma  ashing  can  ba  very  useful  for  the  ranoval  of 
I organic  substrates  and  particles.  The  particles  formed  upon 

* ashing  are  certainly  less  than  1 pm  as  determined  microscopic- 

ally  and  are  probably  considerably  smaller  than  0.1  ym.  There- 
’**  fore,  the  particles  should  be  easily  separated  from  the  larger 

inorganic  particles  originalJy  present.  Oxygen  and  oxygen/nitro- 
4 gen  plasmas  do  extensive  oxidation  to  U02f  probably  forming 

UO^ • In  our  experiments,  about  18%  of  the  UO2  was  oxidized 
to  UO^.  Consequently,  UO,.  does  not  remain  invarient  in  oxygen 

(it  3 ^ 

plasmas . 
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SECTION  IV 


SYSTEM  DESIGN  AND  FLUIDICS 


4.0  Proposed  Schematic  Instrument  Design 

The  schemcitic  instrxament  design  for  a fluidic  particle 
sorter  et  the  time  of  termination  of  breadboard  design  is  shown 
in  Figure  4-1.  Referring  to  this  figure,  the  steps  of  the  pro- 
posed operation  would  have  been  as  follov/s: 

(a)  The  sample  is  placed  in  the  instrument  just  below 
an  aluminum  screen.  The  sample  on  a very  low  ash 
filter  paper  will  be  in  a monolayer  with  particles 
whose  diameters  are  greater  than  5 micrometers 
removed . 

(b)  The.  saT.-.nle'c;  or.i  r-nt.iit.ion  v.ith  respect  to  the  gas 
flov.’  v.’.ill  be  sue))  thcit  the  fi3l  c;r  paper  v/iil  ash 
slowly  from  one  end  to  tl;e  other  releasing  the 
particles  slov;ly  into  tlio  re:  "xining  stages  of  the 
instruinent . Tlie  best  orirntaLion  will  prtbably  he 
that  where  the  norinal  to  . }ie  plane  of  the  filter  is 
perpendicular  to  tlie  gas  flov.’  stream. 

(c)  A miMture  of  15%,  by  volui.^o,  o;.ygcn  in  helium  will 
bo  used  as  the  gas  mixture.  The  total  pressure  will 
be  1 to  2 torr  in  the  ashor  section.  It  should  take 
3 to  6 hours  to  completely  ash  the  filter  paper  and 
introduce  the  particles  into  the  gas  stream.  During 
the  first  30  minutes,  almost  all  the  organic  particles 
present  in  the  sample  will  have  been  oxidized  to 
volatile  products.  Thus,  only  the  interesting 
particles  and  other  inorganics  will  remain  to  be 
carried  off  by  the  gas  stream  to  the  remaining 
stages  of  the  instrument. 
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(d)  The  gas  borne  particles  arc  carried  by  the  gas 
stream  to  the  funnel  design  at  cross  section  A 
v?here  they  arc  collimated  into  sixteen  beams  of 
" ■ 1 mm  in  diameter.  The  helium  flow  at  Gas  Inlet 

B wdll  be  such  as  to  keep  the  flov/  laminar  and 
. thereby  maintain  particle  control.  .The  diameter 

of  the  dielectric  stage  v;ill  be  made  large  enough 
\ to  maintain  a pressiire  differential  of  10:]..  In 

the  ashcr  stage,  the  pressure  will  be  about  1 torr 
I giving  an  oxygen  partial,  pressure  150  microns  of 

Hg,  which  is  sufficient  for  the  ashing  rates  needed 
, to  ash  the  filter  paper  in  3 to  6 hours.  The  di- 


electric stags  v.’il],  have  a diameter  sufficient  to 


maintain  a prer.su.ro  0.1  torr. 

In  the  gas  len.r  created  by  funnel  B,  each  particle 
beam  w.i.1.1.  be  reduced  from  a diameter  of  1000  micro- 
meters to  a ciiamctc.r  of  10  mie:co;'.icterr- . The  gas 
flovj  rates  will  be  adjusted  to  maintain  particle 
control. 

In  tlic  r.pprojcimatel.y  400  micrometer  in  diameter 
and  2.5  to  5 cm  long  cylindrical  consl'.ructi.on,  di- 
electric forces  due  to  the  voltage  on  a 100  micro- 
meter in  d.iamcter  canter  electrode  will  deflect  the 
particlc-;s.  Particle.-'  of  higli  dirl.ectr.l  c constant 
to  density  ratio  will  be  deflected  more  than  those 
of  lov;er  ratios.  Thus,  a ratio  spectrum  is  formed. 

The  microscopic  beam  is  then  magnified  by  the  reverse 
funnel  at  C. 

The  separater  tube  will  select  that  part  of  the  gas 
borne  particle  beam  containing  the  interesting 
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particles,  which  then  pass  into  the  final  stage  of 
separation. 


The  particle  stream  then  passes  through  funnel  D 
where  it  is  collimated  into  a 1 mm  beam.  It  is 
then  diluted  by  gas  entering  at  inlet  C and,  again, 
collimated  by  funnel  E to  a 10  micrometer  beam. 

The  operational  parameters  v;ill  be  found  by  adjusting 
flow  rates  and  diameters  so  that  the  particle  beam 
will  have  one  particle  at  a time  moving  at  constant 
velocity. 


A laser  beam  set  perpendicular  to  the  beam  v;ill 
excite  a fluorescence  from  each  particle  as  it  passes 
by.  An  array  of  filtered  detectors,  each  directed 
at  tlie  intercept  of  the  particle  beam  and  the  laser, 
v;ill  sense  the  f luorescencc . The  output  of  the 
detectors  v.’ill  )‘>e  sent  into  the  control  electronics 
and  pi.ocossed. 


When  there  is  a characteristic  omission  from  an 
interesting  pax'ticle,  the  electronic  processor, 
after  a delay  period,  send  a high  voltage  pulse  to 
a pair  of  deflection  electrodes.  This  delay  period 
allows  for  the  time  it  takes  tJje  selected  particle 
to  travel  from  the  field  of  the  detectors  to  that  of 
the  electrodes.  The  high  voltage  pulse  creates  a 
temporally  short  ion  wind  between  the  electrodes, 
which  deflects  the  chosen  particle  into  a separate 
beam.  If  the  pulse  width  is  60  microseconds,  a 
single  particle  can  be  selected  from  10  particles 
passing  by  per  second. 


(1)  The  deflected  particles  are  collected  on  a filter. 
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4 . 1 The  Original  Design  of  a Precision  Dielectric:  Separator 


With  the  favorable  calculations  of  separation  by  the  cylin- 
drical geometry  for  the  system  of  UO2  and  Si02 , design  of  a prac- 
tical dielectric  separator  was  begun.  This  design  required  the 
location  of  a stream  of  particles  150  ym  from  a 100  ym  wire. 

This  stream  had  to  be  located  initially  at  150  ym  to  ± 5 ym. 

We  believed  that  this  sort  of  precision  and  accuracy  is  possi- 
ble by  taking  advantage  of  streamlines  in  a laminar  gas  flow. 
Figure  4-2  show’s  the  geometry  of  such  a system. 

The  basic  theory  behind  it  is  that  the  particles  of 
interest  are  inserted  In  a large  diameter  tube,  A,  at  a distance 
B betw'een  the  center  line  of  the  tube  and  the  center  line  of  a 
capillary  insertion  device  v/hose  diameter  is  C.  When  the 
particle’  enters  the  gas  flow,  it  introduces  a radial  disturbance 
to  the  parabolic  flov;  pattern  in  the  large  tube.  For  the 
purposes  of  the  follov/ing  analysis  it  is  assumed  that  these 
particles  then  travel  along  the  large  tube  until  the  distance 
L is  reached.  At  this  point,  any  disturbances  to  the  para- 
bolic flow  pattern  introduced  by  the  capillary  tube  is  damped 
out,  and  v/e  again  have  steady  laminar  flov/,  with  the  particles 
having  motion  only  along  the  stream  lines.  This  distance,  L, 
is  probably  on  the  order  of  ten  tube  diameters  or  lOA. 

At  this  point  the  gas  flow'  enters  a converging  nozzle  v/here 
the  flow  is  focused  into  the  dielectric  sorter.  The  dielectric 
sorter  length,  £,  is  sufficient  for  particle  separation  at  the 
chosen  flow  rate.  After  deflection,  the  resulting  particle 
spread  is  magnified  by  a diverging  nozzle  resulting  in  radial 
separation  Y. 

The  first  design  criterion  is  that  laminar  flow  is  always 
maintained  at  critical  points  in  the  system.^  In  the  area  of 


Principles  of  Heat  Transfer,  International  Textbook 
Company,  Scranton,  1965. 
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the  electrodes,  the  geometry  is  at  a minimum.  At  this  point, 
we  must  have  laminar  flow,  therefore  the  Reynolds  Num.ber,  Re, 
cannot.be  over  2000. 


Re  = 2000  = — 
V 


(4-1) 


v;here 


V - mean  fluid  velocity  at  this  point 

a = local  diameter 

V = kinematic  viscosity 

For  this  system,  it  is  assumed  at  present  that  the  carrier  gas 
is  helium,  therefore 

_2 

V = 3 X 10  cm/sec 

By  choofjing  an  appropriate  value  for  "a",  this  yields 
tlje  me.-n  fluid  velocity  at  the  constriction  and  having  our  di- 
electric deflection  calculations  of  the  electrode,  f.  It  also 
defines  the  gas  flow  rate  for  the  system. 

In  this  envisioned  geometry,  the  particles  are  initially 
located  150  -hn  from  the  central  wire,  to  a precision  of  ± 5 ym. 
Considering  that  tlie  particles  v/ill  spread  radially  on  the  or- 
der of  ± 2C  when  initially  inserted  in  the  large  tube  of  diam- 
eter, A,  this  sets  up  a geometrical  criterion  between  C and  B, 
the  distance  between  the  center  lines  of  the  capillary  and  the 
large  tube.  From  this  it  can  be  shown  that  B = 60C. 

We  also  have  a criterion  between  B and  A which  is 


A = 


B 


150 

where  a is  in  microns. 


(4-2) 


I I'ii  n 1^.1,  ■'  Ln.  'PU'^ 


■••wesiw/7.;'  ■^9*fgasT 


In  order  to  get  a feel  for  some  real  numbers,  let  a = 

500  pm  and  C = 200  pm.  These  are  realistic  numbers  con- 
sidering that  capillaries  are  routinely  made  with  I.D.'s  of 
that  size,  and  that  is  fairly  simple  to  machine  a 500  pm  hole 
in  a thin  piece  of  aluminum. 


Using  these  values,  the  velocity  of  the  gas  at  the  elec- 

4 

trode  will  be  3.7  x 10  cm/sec. 


If  one  uses  helium  as  the  carrier  gas,  this  velocity  is 
about  a factor  of  2.5  lower  than  the  velocity  of  sound  in 
helium.  Therefore,  the  system  is  not  subject  to  a normal 
shock  at  the  constriction.  This  velocity  defir  ^s  the  flow 


rate,  v.’hich  v;ill  be  74  cm'^/sec. 


Next  using  the  vnliics  of  C ~ 200  cm,  B ~ 1.22  cm  and 
7i  •-  4.0  cm,  they  define,  a mean  fluid  velocity  at  this  point 
of  5.7  cm,/scc. 


The  Keynolds  nvmibcr  in  tlic  largo  t\ibc  v/ill  be  25. 


This  analysis  proposed  to  sho'w  the  definite  viability 
of  a clielect.rJc  sorter  using  a gas  flow  to  locate  the  par- 
ticles. There  are,  however,  design  problems  that  must  be 
looked  into  and  experiments  that  would  have  had  to  be  done  in 
order  to  develop  the  optimal  design  parameters. 


With  the  sample  numbers  used  in  this  analysis,  the 


maximum  allowable  fluid  velocity  at  i.hc  electrode  comes  out 

4 

to  be  3.7  X 10  cm/sec.  This  high  velocity  is  undesirable 


for  several  reasons.  One  is  that  compressible  effects  have 
to  be  taken  into  account  at  such  high  velocities;  two,  since 
the  sorting  deflection  takes  place  on  the  order  of  milli- 
second this  implies  an  electrode  length  of  about  one  foot. 
This  is  impractical  with  respect  to  the  hole  size  that  would 
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have  to  be  drilled  into  the  piece  of  metal.  One  could  lower 
the  flow  rate  and  thus  the  velocity  at  the  constriction,  let 
us  say  a factor  of  ten  or  twenty » ie. 


V = 1830  cm/sec. 


This  would  give  a reasonable  velocity  at  the  electrode  and  a 
reasonable  pressure  at  the  electrode  (assuming  one  starts 
with  approximately  100  torr  in  the  large  tube) , but  the  flow 
rate  would  now  be  only  3.7  cm^/sec.  The  Reynolds  number  in 
the  larger  tube  would  only  be  1.25.  At  tl.is  low  Reynolds 
number  an  experiment  would  Jmvo  to  be  done  to  see  if  there 
is  still  a parabolic  velocity  distribution  in  the  large 
pipe.  At  these  lov/  velocities,  temperature  control  might 
be  a proJjlem.  An  am  lysis  v;ould  have  to  be  done  on  the 


convector  liei’ti)ig  effects 

Tlu.‘  ext,ct  end  effects 
ViouDd  1„  vr  to  be  ei:.^.i  rice.' 
to  go  into  the  optimum  v;ay 


along  the  tube  length. 

Oi  t]je  Cc'i'ij.loiy  ill  tj.i.:  flow 
.ly  asi'iortaii’f  d . A study  \,ou]d 
of  inserting  the  partieles. 


field 

linve 


Various  pliysir-el  propc-rt  ic.::.  of  the  system 
be  experimentally  measured  a)iu  dealt  with,  i.c. 
mode-  of  the  central  wire  or  tJ.e  acoustic  resona 


\;ould  liavo  to 
the  vibration 
nces  of  the 


large  tube. 


This  flow  and  particle  placement  problem  is  a complex  one, 
but  it  v;ould  be  soluble  if  the  dielectric  sorter  did  not  require 
such  low  total  pressures.  The  experimental  mock-up  as  outlined 
could  have  been  readily  adapted  for  prototype  use.  For  example, 
in  a prototype  sorter,  many  more  than  one  capillary  could  be 
introduced  radially  in  the  gas  flow.  Perhaps  a continuous, 
concentric  tube  device  could  be  used  instead  of  capillaries. 
Thus,  significantly  decreasing  the  time  to  sort  all  of  the  par- 
ticles in  question. 
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The  most  difficult  construction  problem  would  have  been 
the  exact  alignment  of  an  approximately  12  cm  long,  100  um  in 
diameter  wire  inside  a 2 to  5 cm  and  500  m in  diameter  tube. 
Since  the  alignment  accuracy  must  have  been  to  ± 3 um  while  hav- 
ing very  high  tension.  Very  high  tension  is  required  to  reduce 
sensitivity  to  vibration  and  other  deformations. 

4.2  Pro£o^s^P2six;n_of  Determine  Minimum  Operatinq 

Pressure  for  Particle  Control 

In  an  effort  to  ascertain  an  optiiuum  design  parameters  for 
the  fluid  sorting  system,  the  eeuatrons  of  Section  4.3  that  re- 
late the  various  para^noters  in  Figure  4-2  are  plotted.  Figure 
4-3  was  generated  by  relating  the  size  of  the  constriction  to 
the  other  parameters  such  as  the  large  tube  diameter.  A,  con- 
sisting of  a fainily  of  curves  depending  on  the  insertion  tube 
diameter.  C;  the  toLd  flow  r.-.to,  Q,  th^nugr.  the  system;  the 

fiuid  velocity  at  tl.e  conrdriction,  v,  and  tie  constriction 
lengtli,  3. 

^An  orpr.r>ment  was  proponed  by  the  use  of  these  graphs  as 
a rercrenco.  A constriction  size  diameter  of  500  um  was  chusen. 
This  dictated  the  fluid  velocity  at  the  constriction  to  be 
1/07  c;m/scc.,  the  constr id  ion  length  to  be  1.8  cm  and  tlic  flow 
rate  to  be  2 1C  cc/min.  Those  seemed  all  to  be  reasonable  dimen- 
sions to  worh  with.  Since  the  asher  was  tentatively  designed  to 
avc  a di.mirtcr  of  7.C  cm,  this  was  chosen  for  the  fluid  control 
experiment  also,  dictating  an  insertion  tube  inner  diameter  of 
381  um.  This  also  constricts  the  position,  B,  of  this  tube  to 
be  2.3  cm  from  the  centerline  of  the  largo  tube  A. 

The  next  question  to  have  been  addressed,  was  what  was  the 
local  velocity  of  the  carrier  fluid  at  point  B.  Using  a para- 
bolic fluid  distribution,  this  was  found  to  be  0.01  cm/sec 

Ideally  one  would  want  to  match  this  to  the  insertion  velocity 
of  the  particles. 
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ft  ELECTRODE  LENGTH-C:/. 

1;  = FLUID  VELOCITY  AT  ELECTRODE  Cf/i/SEC 
Q = FLUID  FLOW  RATE  CC/SEG 
A = LARGE  TUBE  DIA-CI/i 

a=  electrode  dia-cm 

C = INSERTION  TUBE  DIA-CM 


FIGURE  4-3 

PLOTTING  OF  FLUDIC  PAR7VMETER 
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Next,  a calculation  was  done  to  determine  the  terminal 
velocities  of  uranium  oxide  particles  and  silicon  dioxide 
particles,  varying  from  5 to  0.1  -,jm,  if  the  forces  acting 
upon  them  were  just  gravity  and  a fluidic  drag.  It  was 
found  that  the  terminal  velocities  varied  from  6.4  cm/sec 
to  0.058  cm/sec,  where  ideally  one  is  looking  for  insertion 
velocities  in  the  neighborhood  of  0.01  cm/sec  nominal.  This 
means  that  first  of  all  any  expcrimo}it  v;ill  have  to  be  posi- 
tioned vertically,  so  that  the  insertion  velocities  v.’ill  be 
minimised.  II  also  means  that  one  will  not  be  able  to  match 
particle  insertion  velocities  to  Isulk  fluid  velocities. 
Hopefully,  ol  ' tills  v.’ill  mean  is  that  the  particles  could 
separate,  but  along  one  streaiulino.  If  mounted  vertically, 
this  should  be  an  unimportant  effect. 

Anoth''-"  calculatic.il  v.ns  perjormod  in  to  .atimatc 

the  perturlintion  effect  of  the  capillaiy  insertion  tube.  It 

Wins  found  th.il  the  pcj  erjit  }>e.vtu”b;ition  to  tlic  flov.’  field^ 

2 

is  a function  of  (^q/L)  ' where  9.^  is  the  length  of  the  in- 
sertion tubo  .-’nci  L is  a length  '.uasured  frcni  the  tube, 
dov.-nstrean.  For  exru;iplc,  if  one  v.’ants  a one  percf  it,  or 
less  pertu  :bi  ' ion  fron  tlto  capillary  and  tli''  capillary  is, 
let  us  say  10  cm  long  v.’hirh  yields  a length,  L.  of  100  cm. 
This  was  to  )io  used  as  a design  critcric^n  for  the  length  of 
the  larger  tube  of  diameter.  A,  before  the  funneling  begins, 
and  thus  the  focussing  of  the  fluid  streamlines. 


1.  "Theoretical  Hydrodynamics",  by  Mi line -Thoms on, 
MacMillan,  1961. 
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• 3 The  Effects  of  Brownian  Motion  on  Particle  Positioning 


When  the  particle  diameter,  d,  is  much  greater  than 
the  mean  free  path  of  the  gas,  X, 


X = (2kTt/37idn) 


(4-3) 


where 


X is  the  mean  displacement  perpendicular  to  the  direction 
of  flow  of  a partir-le  in  a beam  during  a time  t, 

k is  the  Eoltzman  constant, 

T is  absoJ.ute  temperature,  and 

ri  is  the  viscosity  of  tlie  gas. 

Then  the  beam  diameter,  D,  after  a time  t is  given  by 


D = H-  2X 


(4-4) 


For  ]k:  £.  i-  JOO  torr,  hov.’evc.r , the  mc*aM  free  pc-.tl)  is  1.3  yrn. 
Since  this  is  not  small  compared  to  the  particle  sizes  (5  ym 
to  0.1  tim)  , the  empirical  formula  due  Millikan^'  can  be  used 
to  get  t)jc  effective  value  of  n by  divn cling  ]:>y  the  factor 


1 1 (2X/d.)  (A  -I-  E exp  [-Ccl/kX]) 


(4-5) 


\’herc 


A = 1.23 
B = 0.41 
C = 0.88 


(4-6) 


For  a particle  beam  moving  at  5.73  cm/sec ^ down  a tube  of 
length  38  cm,  the  travel  time  is  6.63  sec.  The  number  was 
obtained  from  our  preliminary  design  in  the  previous  sections 


1.  Millikan,  Phys.  Rev.  1 (1923) 

2.  This  value  should  be  taken  as  an  upper  limit  to  the  velocity. 
It  has  been  determined  that  slower  rates  would  be  more  desir- 
able from  a fluidic  point  of  view;  a slov;er  rate  would  give 
an  even  greater  beam  spread. 
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and  represents  the  type  of  condition  that  would  be  present  in 
an  automatic  fluidic  sorter.  To  determine  X,  the  mean  free 
path,  the  following  relationship  is  used. 

X = kT/(2^7rc^P)  (4_7) 

where  c is  the  molecular  (atomic  dimeter)  and  P is  the 

O 

pressure.  For  helium,  c = 2.20  A.  For  an  initial  beam 
diameter  of  1.0  mm,  the  calculated  beam  diameter  at  the 
end  of  tube  is  shown  in  Table  4-1  for  different  particle 
sizes  and  at  various  pressures.  Since  the  particle  distri- 
bution is  not  uniform  and  probably  has  a Gaussian  distribu- 
tion, the  diameter  represents  that  cross-section  where  95% 
of  the  original  particles  will  pass  through.  These  results 
show  that,  at  pressures  greater  than  10  torr,  there  is  rela- 
tively good  particle  control.  Hov.’ever,  when  the  beam  diver- 
gence is  greater  than  50%  (i.e.  1 mm  to  1.50  mm),  there  is 
insufficient  particle  control  for  proper  operation  of  dielec- 
tric sorting  and  one-at-a-time  sorting. 

TABLE  4-1 


THEORETICAL 

SPREAD 

IN  A PARTICLE  BEAM  DIA.METER 

P 

X (^im) 

Particle  1 
10  pm  5 pm 

Diameter  (mm) 
1 pm 

0.5  pm 

0.1  pm 

4 atm 

0.047 

1.01 

1.02 

1.04 

1.06 

1.17 

1 atm 

0.19 

1.01 

1.02 

1.04 

1.07 

1.29 

100  torr 

1.44 

1.01 

1.02 

1.08 

1.16 

1.77 

10  torr 

14.4 

1.03 

1.05 

1.25 

1.49 

3.44 

1 torr 

144 

1.09 

1.15 

1.77 

2.54 

8.72 

0.1  torr 

1440 

1.30 

1.49 

3.44 

5.88 

25.42 

4.4  Conclusions 


In  section  2.3,  a pressure  of  0.1  torr  was  determined  to 
be  the  highest  pressure  tolerable  for  good  dielectric  separation. 
However,  in  this  section,  it  was  theoretically  determined  that 
there  is  insufficient  particle  control  for  particles  1 ym  or 
less.  Consequently,  for  particles  1 ym  or  less,  dielectric 
separation  is  not  possible.  In  addition,  during  ashing  of  the 
substrate  to  make  the  particles  gas  borne,  high  gas  pressures 
and  velocities  are  required  to  prevent  the  particles  from  stick- 
ing to  the  sample  holder.  Although  the  exact  pressures  and 
velocities  were  not  determined  by  experiment,  they  T/ill  be  at 
least  two  orders  of  magnitude  greater  than  the  maximums  calcu- 
lated for  dielectric  separation.  The  theoretical  irreconcil- 
able differences  in  our  automatic  particle  sorter  design,  which 

IS  the  only  possible  design  using  the  techniques  considered, 
are  sumniarisod! 

(1)  nigh  pressure  and  high  flow  rate  are  required  to  get 
the  particles  into  a gas  borne  state.  Preliminary  experimental 
investigation  showed  that  this  transfer  from  substrate  to  gas 
flow  is  not  likely  to  be  with  low  loss. 

(2)  High  pressure  and  lov;  flow  rate  are  required  to  posi- 
tion the  particles  to  ± 5 ym  in  the  dielectric  sorter  and  one- 
at— a-time  modules. 

(3)  I,ov/  pressure  is  required  to  overcome  particle  size 
dependence  and  produce  a sufficient  dielectric  constant  to  dens- 
ity ratio  spectrum  for  a bulk  separation. 
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APPENDIX  A 


BASIC  DIELECTROPHORETIC  SEP  iPATIQN. 

theoretical  developme n t 

A.  0 Introduction 

It  IS  well  knov.'n  that  an  uncharged,  dielectric  particJe 
placed  in  a nonuniforn  electric  field  experiences  a net  force 
tending  to  move  the  particle  towards  regions  of  higher  field 
intensity.  In  this  appendix,  this  effect  is  analyzed  for 

Its  application  to  the  separation  of  particles  in  the  0.1 
to  5 ym  range. 

The  dielectric  force  on  a particle  is  evaluated  as  a 
function  of  its  dielectric  constant,  density,  eto.  and  the 
applied  electric  field,  for  tvio  represontativo  electrode  con- 
figurations, and  the  results  for  the  two  cases  are  compared. 
Later  it  IS  shoun  that,  in  some  sense,  these  configurations  are 
optimally  effective  in  prodiK.inn  dielcetrophoi-esis . 

Nett,  tlic  effect  of  the  pre.senoo  of  a fluid  medium  on  the 
eotual  movement  of  partielos  is  analysed.  The  viscous  forces 
on  particles  of  various  size  are  evaluated,  due  account  being 
taken  of  their  extrciely  small  diameter.  The  actual  movement 

o a particle  is  then  derived  and  described  in  terms  of  a nor- 
nialized  equot.lon  of  motion. 

Finally,  a number  of  criteria  are  derived  vhioh  are 
effective  in  determining  the  efficiency  of  dielectric  separa- 
tion. These  critoria  are  of  use  in  the  actual  design  of  a 
pairticlc  separation  stage. 
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3.1  Dielectric  Acceleration  of  Particles 


I, 


The  surface  density  of  charge  induced  on  a dielectric 
particle  (of  dielectric  constant  K)  placed  in  an  electric 
field  of  magnitude  E (refer  to  Figure  A-1  for  the  geometry)  is 

(A-1) 

where  eo  is  the  permeability  of  free  space.  This  result 
follov/s  from  Gauss'  Law,  relating  the  electric  displacement 
to  induced  surface  charge. 

Assuming  a irod  like  shape,  A is  the  cross-sectional  area 
of  the  particle.  The  total  charge  on  one  face  of  the  particle 
is 

q - AE(,  (K  ~ 1)E.  (A-2) 

With  a pa.vticle  lenglh  %,  tlic  moment  of  t)ie  induced  dipo].e 
is  tlien 


1 i 


\ 


\ 


q£  = AEpf  (K  - 1)E  = £5V(K  - 1)  , 


(A-3) 


where  V - h%  is  tlie  particle  volume. 

Althougli  this  derivation  assumes  a cylindrical  particle, 
the  dipole  moment  and  other  quantities  to  be  derived  are  inde- 
pendent of  the  actual  particle  shape  and  its  spatial  orientation. 
h more  general  analysis,  using  the  relation  that  the  induced 
surface  charge  is  proportional  to  the  normal  component  of 
field  and  then  taking  a surface  charge  integral,  proves  this 
independence. 

The  following  analysis  thus  applies  to  particles  of 
arbitrary  size,  shape  and  orientation.  We  will  see  that  the 
acceleration  of  a particle  is  independent  of  those  parameters. 

' 5 

I 
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Let  the  radial  dependence  of  the  electric  field  be 
expressed  as  L (r)  along  the  normal  to  the  plane  electrode  of 
Figure  A-1.  For  a particle  at  radius  r,  the  forces  on  each 
face  due  to  the  induced  cliarges  are  as  follows: 

Force  on  tlie  negatively  charged  face  is 


F_  = -a^AE(r)  . 


(A-4) 


Force  on  the  positively  charged  face  is 

F^.  - a.AE(r  •!-  £)  . (A-5) 

T1j.cj  net  force  on  tlje  particle  is  therefore 
E r=  E_  ,,  a.A{E{r  -i-t)  - E(r)  }. 


In  a uniform  field,  i.e.,  AE  = 0,  F - 0.  Thus 
net  force  only  in  nortuniform  fields,  i.e.  AE 
E(r  •!-  1)  may  lie  oxprers.ed  as  a function  of  E(r) 
a Taylor's  serie.s  expansion  aboul.  r: 


there  is  a 
0 . N ov/ , 
by  means  of 


E(r  + f)  - E(r) 


t?E/9r  + H. 

Dr^ 


• (A- 7) 


It  is  assumed  that  the  second  and  higher  dcj.'ivatives  of  E 

can  be  neglected  with  respect  to  9E/Dr.  Or,  more  specifically, 
it  is  assumed  that 

<<  2 ^ ^ 

9r^  I 8r  (A-8) 
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This  assumption  is  justified  since  the  particle  sizes 
are  small. 


/,  E(r  + H)  = E(r)  + £3E/Sr. 


(A-9) 


The  net  force  on  the  particle  is 


F = a, A f-4-—  = e (K  - 1)  VE  , 

o 9r  ' 


3r 


(A-10) 


by  substituting  Equation  (1)  into  (6) . 

In  free  space ^ the  acceleration  of  such  a particle 
(density  p)  is  given  by 


a 


field 


_F_ 

m 


PV 


So  (K 


1)  VE 


3E 

9r 


o 


(K  - 1)  E 


3E 

3r 


(A-11) 


To  determine  the  field  E between  the  charged  wire  and 
the  plav.e  electrode  of  Figure  A-1,  the  method  of  image  charges 
is  used.  VJe  assume  the  presence  of  an  image  line  charge  of 
opposite  sign  an  equal  distance  below  the  plate.  The  surface 
of  tl-ie  plute  is  then  equipotential , and  the  two  configura- 
tions give  rise  to  identical  fields.  it  is  therefore  suffi- 
cient to  derive  the  field  expression  for  a pair  of  wires  of 
radius  ro  separated  by  a distance  2r^. 

By  Gauss'  Theorem,  the  field  intensity  E^(r)  produced 
by  the  upper  wire  is  related  to  the  charge  on  it  by 
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J 


Q = + = + egE^(r)  2Tir  (A-12) 

or, 

E (r)  = — . (A-13) 

2-nc  r 
o 

Similarly,  tho  field  component  due  to  the  lov^er  vrijre,  E (r)  , 
can  be  shov.ni  to  be 


E (r)  - 


2tte 


Q 

Ur'^  -'iT  * 


(A-14) 


Ujjing  the  su) 'or]r»osition  p?'inciple. 


the  total  field  is  then 


E(r)  •---  E , (r)  -h  E_(r) 


(A-15) 


Equation  15  cdui  nov;  be  solved  in  terms  of  the  total  voltage 
botv.’een  the  upper  voire  and  the  g3.ound  plane.  Let  this  voltage 
be  Vg . Then 


V 

0 


(A-16) 


^’he  integrand  of  Equation  16  can  be  evaluated  to  obtain 
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,.-U.  - : 


dr  1 , 1 

r - 2r^r  ^ 


1 ^ 


Therefore , 


Q r^i  “ ^ 

V = ^ — In  — ^ 

0 2tj€^  \ r 


(A-17) 


(A-18) 


This  equation  relates  the  charge  per  unit  length  on  the  v/ire 
conductor  to  its  voltage.  Thus  the  capacitance  per  unit 
length  of  this  geometry  is 


^ {2r,  *-  r ) 

C - 5-  2 Ko/hi  — I “ 

Vo  r 0 


(A-19) 


Let  In  (2r^^  r /r  ^ be  defined  as  P,  a dimensioailess  con- 
stant. Then,  from  eque-.tion  18, 


_ OP 


(A-20) 


27re^Vo 


(A-21) 


We  substitute  (21)  into  equation  (15)  to  obtain 


E(r) 


(2r^  - r)  r ' 


(A-22) 


From  the  above  equation  it  is  possible  to  obtain  3E/9r  as 
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8E(r) 

3r 


3 

3r 


(2r^r  - r^)~^j 


ill  (2im-  - )~^(2r, 


2V.r 

1 (21, 

P 1 

4V^r^  (r  - r^) 


- 2r) 


(2r^r-i-2)^ 


and 


3^E(r)  ..  1 

-2  F 2 


IbV^ri 


(r,  - 


r^-) 


Dr 


(21-3^1-  ~ r'^)  2 


(2r^r 


2..  2 


E(r)  _ 


.•.  E r rr  .. 

3r  p 

P" 


(r  - r^) 


(2r.r  --  r ) 


2,  3 


ThuG,  from  Equation  (n) , the  net  aoceleration  of  the  par- 
tide  in  free  space  v.’ill  ]:>e 


8e 


"field 


-i  (K  - 1,  FF!  . Jld  t 

^->v-  „2,  3 


(2r  r - r ) 


A.  3 Detorijilnatio.:  ^ ||  for  a Cyli  ndrical  Geotletrn, 

The  dielectrophorotic  effect  is  calculated  hero  for  the 
alternative  cylindrical  geometry  shown  in  Figure  A-2.  Only  a 

portion  of  the  outer,  completely  cylindrical  conductor  is 
shown  in  this  figure . 


The  net  force  on  the  particle  is  again  given  by  equation 

( J.0 ) • 
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(A-23) 


(A-24) 
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FIGURE  A-2 


CYLINDRICAL  ELECTRODE  CONFIGURATION 
FOR  PARTICLE  ACCELERATTDNl 


AND  DIMENSIONS 


I 

I 

I 

I 


:! 


DE 


Dig  E tcriii  must  I^g  doterrainou  for  cylindrical 
geometry  by  Gauss'  Theorem  as  sliown  by  equations  (12) 
and  (13) . In  this  case,  there  is  only  a single  field  component. 


E (r)  = E,  (r)  = ^ — . 
+ 2nc,r 


(A-27) 


Equation  27  is  .solv(-d  to  oliLain  tlio  total  voltage  bet\7een 
the  outer  and  inner  conductors.  Integrating  Equation  2 7, 


V. 


E(r)  dr 


•i ^ Q 


2ire 


In  r 


(A-28) 


Thus,  in  this  cnase,  the  cay'acitance  C per  unit  length  ii 
C = 


_ 0 .. 


'o  ln(r^/ry)  • 


(A-2i') 


Let  lji(r^/rj  be  defined  nr,  P , a dimension Ic.s.s  con.stent. 
Thc:n , } rom  this  equation, 


V, 


QP 


2TTt 


(A- 30) 


or 


27re 


Q = 


0 

r-  V, 


(A-31) 


We  substitute  Equation  (31)  into  Equation  (27)  to  obtain 


V. 


E(r)  = 4 I . 

P ^ 


(A-32) 
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By  differentiating  Equation  (32) 


JE  _ 1 


(A-33) 


a^E 

ar2 


2V 


0 1 


P r 


3 ' 


(A-34) 


and  thus 


9E 

af 


'2  3 

P r-* 


(^-35) 


Prom  Equation  11,  the  net  acceleration  of  the  particle 
in  free  space  becomes 


c . (K  - .1 ) V 


Afield 


c 1 


(A-3G) 


A.  4 Com  r:i  r-on  of  j j -'~t  jj-aj  and  Planar  Cer>r,-.r* trier 

On  the  basis  of  the  prcvlour.  derivation,  tlie  tv:o  geo- 
ineti'ics  can  bo  coiTiparf;d  to  dotoimine  their  relative  effi- 
ciencies in  producing  dielectrophoretic  acceleration . The 
two  key  equations  to  bo  coup:  I'ed  are 


“field 


P r.2 


(2r^r  - r^) ^ (A-26) 


and 


Vo‘ 


o field  " ^ “T  . (A-36) 

K P r 

Let  X be  a normalized  position  parameter,  such  that 
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X - r/r^  , 


(A-37) 


where 


1 i 1 . 


(A-38) 


Then  Equation  (26)  and  liquation  (36)  are  normalised  an  fol- 
lov7C : 


Afield 


E,(K  - 1) 

3 ^0  , (A-39) 

pr^  P'^(x  - hx^)^ 


and 


" lield 


Cc(K  - 1)V 


pr, 


0 .1 


(A-40) 


The  field  accel<;rationc  indneed  by  tlic  tv;o  electrode 
co/i f i.gurationo  are  identiccil  functionn  of  the  applied  vol- 
tago,  Ij.near  di:..enG:i.ons  of  the  geonK  luy,  etc.  Tt  in  suffj- 
ciont  to  compare  tlie  ternio  P and  P , and  tlie  factors  involving 

I 

X,  rn  order  to  comg'urc  a vrith  e . 

If  the  two  geometricri  can  be  considered  identical  v/hen 

tlie  same  i.n  both  cfi.ses,  then  v;e  can  compare  P v’ilh 

I 

P . 

(2r,  - ro) 

, P = In  in  {(2r^/r„)  - 1}  (A-41) 
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P = ln{r^/rp) 


(A-42) 


By  definition  is  always  > 1, 

1 ^ (A-43) 

and,  since  .ln(;-:)  is  a positive  iTionoton;i  c function  of  >:,  j.t 
folJov.’s  that 

I 

P > P also.  (A-44) 

On  tlic  basis  of  coi.iparison  of  P and  P'  alone,  the  accel- 
erations produced  i^y  flu;  cylindrical  ueuir.etry  will  be  larger 
than  tlio  accelerations  produced  by  the  planar  configuration 
of  clcctrcdcs. 

Tiie  otlicr  basis  on  v;hich  tlieso  accel-orations  r.,ay  be 
coi’iparcd  xs  tlic  ratio  of  the  functions  involving  x In  the 
equations  J9  and  40.  This  is  a nunerlcal  coraparison , and  is 
carried  out  in  tlic  next  seclion.  I\  r the  present  purposes 
it  i suificiont  to  state  that  tlrr  factors  involving  x 
diffej  only  .sliglitjy  from  one  another,  and  thus  on  this  basis 
the  two  geometries  arc  almost  identical. 
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A . 4 . 1 Compari.son  of  the  Cylindrical  and  Planar  Geometries 
Based  on  the  Factors  (x)  and  (x) 

It  was  determined  in  Section  A. 4 that  the  cylindrical 
and  plc'inar  electrode  confiejurations  can  be  compared  on  the 
basis  of  the  functions 


F^(x)  - 


(1  - X) 

( i 3, 

(x  - 2>'  ) 


and 


^2 


1__ 

3 

X 


Then  a measure  of  the  relative  efficiency  of  the 
cylindrical  geometry  over  the  planar  is  the  ratio  F2(x)/i'j(x). 

The  factors  F^  (x)  , F2(x)  and  tho.i  r ratios  v.’ere  computed 
for  the  c.ntire  rarjrjc  (0  to  1)  of  x in  increments  of  0.01.  The 
following  table  (Table  1)  lists  values  of  (x)  , F^Cx),  and 
their  ratio  for  various  values  of  x,  tlie  normalii^ed  position 
parameter . 

It  can  be  concluded  fro;u  Table  A- 1 that  there  is  no  signifi- 
cant cid'/;;,ntage  t.o  be  obtained  from  selcctijuj  one  of  these  geo- 
metries over  the  other. 

The  computer  program  used  to  genei ate  F^  (x) , F2 (x) , and 
their  ratio  is  shown  in  Figure  A-3. 
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Comparison  of  the  Cylindrical  and  Planar 
Geometries  Based  on  the  Functions' F, (x)  and  F^ (x) 
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Figure  A- 3 
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A.  5 Effects  of  the  Carrier  Gas 


So  far  the  viscosity  of  a carrier  stream  gas,  which 
is  necessary  to  move  the  particles,  has  not  been  considered. 
The  behaviour  of  particles  under  dielectrophoretic  forces 
is  considered  in  this  section  in  the  situation  where  the 
particle  motion  occurs  in  a gaseous  medium.  The  medium 
will  be  considered  to  be  at  rest  (i.e.,  to  have  no  radial 
velocity,  although  a transverse  stream  flow  velocity  does 
exist) , and  the  particle  is  moving  through  this  viscous  medium 

The  effect  of  viscosity  on  the  motion  of  very  small 
particles  is  described  by  Stoke 's  Law  (shown  later) . The 
viscous  force  exerted  on  a particle  moving  v;ith  a velocity 
V is  directly  proportional  to  the  particle's  diameter. 

.Since  the  mass  of  the  particle  is  proportional  to  the  cube 
it®  diameter,  the  ratio  of  force  to  mass  for  a given 
velocity  depends  inversely  on  the  square  of  the  particle 
diaiii'itcu  . Thus  for  very  small  di  arn':  •or  particles,  visco- 
sity is  capable  of  inducing  very  high  d^.g  aecclorationp . 
Consequentially,  the  effect  of  viscosity  is  very  important 
in  the  case  where  small  particles  are  being  accelerated 
in  a gaseous  medium. 


A 5 . 1 Stoke 's  T.aw 

^orce  required  for  steady  motion  of  a sphere 
through  a viscous  mediiam  according  to  Stoke 's  Law  is 

P = 3irndv,  (A-45) 

where  d is  the  particle  diameter  and  n is  the  viscosity 
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of  the  medium.  It  is  assumed  that  the  velocity  v of  the 
particle  is  relatively  small.  This  statement  of  Stoke 's 
Law  is  in  its  simplest  form. 


This  formula  is  valid,  hov/ever,  only  for  the  range 
of  values  for  which  the  diameter  of  a falling  sphere  is 
much  larger  than  the  molecular  mean  free  path  of  the  gas 
medium.  Since  a typical  value  of  the  mean  free  path  for 
air  at  STP  is  0.06  ]im,  this  form  is  accurate  only  for  the 
larger  particles  of  the  0.1  to  lOpm  in.  diameter  range. 

The  intermediate  condition  in  v;hich  the  particle  size 

approaches  the  rnolecvlctr  mean  free  path  leiigth  is  not  well 

analyzed  theoretically.  In  this  region  the  empirical  for-- 

★ 

mula  of  Mi.llikan  is  knovm  to  be  almost  as  accurate  as  the 
mac ros  cop i c S t oke  ' s Lav; . 

Millil'.r.n's  empirica]  law  is  as  f(.illov;s: 


r =-•  — 


3ti  ndv 


1 H I ja  -1 


(A-4C) 


v.’hero  is  iilie  molecu.lc:r  mean  l.rcc  pat.li  of  the  viscous  gas , 
and  a,  b and  c are  three  empirically  determined  dimension- 
less ccjista.jits  v;ith  the  following  values: 


a = 2.46, 
b = 0.82, 
c = 0.44. 


(A-47) 


Note  tliat  Millikan '.c  equation  in  the  limiting  case  of  L/d 
<<  1 reduces  to  Stoke 's  Law  for  large  particles.  Milli- 


Millikan,  Phys . Rev.  22  (1923),  1. 
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kan's  equation  for  the  drag  force  on  a particle  is  used 
in  the  subsequent  analysis.  The  advantage  of  using  this 
is  that  the  motion  of  the  particle  can  be  determined  for 
various  pressures  (and  therefore  moan  free  path  lengths) 
of  the  gas  stream.  The  denominator  of  A-46  gives  in  effect 
the  low  pressure  correction  for  which  is  then  over  the 

operational  range  of  gas  pressures  independent  of  the  gas 

* 

pressure. 


A.  5 . 2 Equation  of  Motion 

Since  the  viscous  drag  force  on  the  moving  particle 
is 


^drag  " 


(A-4  8) 


the  deceleration  of  the  particle  due  to  viscous  drag  xs 


a 


drag 


^drag  _ ‘ 

pd^  + pdL(a  + be”^'^''^) 


(A-49) 


For  a particle  of  diameter  d suljjected  to  a dielectrophor- 
etic  force  and  moving  in  a viscous  medium,  the  net  accel- 
eration is  then  given  by 


a 


net 


Afield 


(K-l)Vo 


drag 

2 


,2  3.3 

1 


pP  '“r,  'X 


18n  ar/9t 

2 ~ / V-  -cd/L.‘ 

pd  + pdL(a  + be  ' ) 


(A-50) 


In  this  equation,  the  cylindrical  electrode  configuration 
has  been  assumed. 

* Sec  "Kinetic  Theory  of  Gases,"  Kennard,  McGraw-Hill, 
New  York,  1938. 
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A. 5 . 3 Normalized  Equation  of  Motion 


In  terms  of  the  norimilized  position  coordinate  x, 

I " 

and  tlia  normalized  velocity  x and  accleration  x,  the  pre- 
vious equation  can  be  written  as 


rix 


c^CK-DVg^  18nr^x 

~rT2  T3  ,2  ^ . , -cd/L.  ' 

pP ' 1^2  ^ pd  + pdL(a  t-  be  ) 


(A-51) 


or,  dividing  throughout  by  r^,  we  have  the  normalized 
equation  of  motion  of  the  particle  as 


II 

X 


e,(K-l)V„ 

’ T,  . 2 4 3' 

pP  ' X 


18nx 


pd^  + pdL(a  + be 


(A-52) 


This  equation  together  v’ith  the  state  equations 


and 


I I 

V = r - r^x 


net 


■I 


r 


1 


(A-53) 


(A-54) 


define  the  moVion  of  a dielectric  particle  moving  under  the 
influence  of  a diclcctrophoretic  force  in  a viscous  medium. 


A. 6 Upper  Limits  on  the  Term  E used  in  Computing  the 
Dieioctrophoretic  Force  on  a Particle 

So  far  in  this  analysis,  only  two  out  of  a large  num- 
ber of  possible  geometrical  configuration  of  electrodes 
have  been  considered.  Although  both  of  these  geometries 
were  shown  to  bo  equally  effective  in  producing  a nonuni- 
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form  field  it  remainf*  to  be  shovm  that  more  effective 
configurations  do  not  exist.  Accordingly  at  this  point, 
it  is  preferable  to  talce  a more  general  approach  and  to 
consider  an  argument  which  sets  an  upper  bound  on  the 
variation  of  the  function 

with  r.  It  is  possible  to  evaluate  such  an  upper  bound, 
based  on  fundamental  electrostatic  principles,  that  is 
independent  of  any  particular  electrode  geometry'. 

To  evaluate  this  maximum  one  must  consider  the  total 
electric  field  function  produced  at  a point  r in  space  due 
to  a disci'ete  distribution  of  charges.  In  a physical  situ- 
ation, the  distribution  of  these  charges  will  be  continuous; 
however,  no  loss  of  generality  is  incurred  in  considering 
a discrete  charge  di sc’ibution . 

It  is  convenif.-nt  to  consider  two  possible  classes  of 
electrode,  configurations  separately;  one,  in  v;hich  an 
invariance  of  the  electrode  geometry  exists  along  one  o.'^  the 
dimensions.  This  is  the  case  with  the  cylindrical  and  planar 
geometries  tliat  have  been  considered.  The  second,  and  more 
general  ca.‘'e,  is  one  in  which  the  electrode  configuration 
varic.s  arbitrarily  along  all  thiee  dimensions.  The  elements 
of  tlic  first  configuration  are  units  of  line  charge,  of  the 
second,  point  charge  units.  V\’e  will  call  the  first  configu- 
ration the  two  dimensional,  and  consider  this  first. 


A. 6.1  Two-dimensional  Geometries 

If  a number  of  line  charges  located 
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at  points  r , r„  r then  the  net  electric  field  at  a 
1 -c'  ' n 

•point  r in  space  due  to  these  lino  cliargos  is  the  vector 
sum  of  the  field  duo  to  each  such  charge. 


E(r)  = 


1 I 


'n 


Avc 


(A-55) 


from  this , 


3r 


_1__ 
A It  c 


+ 


(r 


(r  - i-j) 


-I t 


‘n 


(r  - r^) 


:a-56) 


7W  9E 

Now,  both  E(r)  and  ~ and  therefore  E have  significant 


va.Tu';  only  at  tlic  sing\:.1  a/ii  les 


r = ^2,  • ••  or  r^. 


Also,  at  each  of  those  po.-intr.,  the  contribution  due  to  all 
other  line  cliorgcs  is  relatively  insignificant.  This  can 
be  seen  from  the  fact  that  both 


'n 


Ti~+  c - r J 


(r^  + e rp"  ’ ^ 


and 


(r  + c - r,) 
1 1 
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Therefore  both  the  field  and  the  gradient  of  the  f;.eld  are 
a maximum  in  the  vicinity  of  a line  charge.  It  is  thus 
sufficient  to  consider  the  maximum  value  of  E in  the 

0 jT 

vicinity  of  a single  line  charge  (for  convenience,  at  r - 0) . 
The  field  due  to  this  charge  is  expressible  as 


E(r)  =3. 

r 


(A-57) 


and 


(A-58) 

(A-59) 


At  this  point  v/e  can  compare  our  result  with  Eqr. . 25 

and  more  directly  v;ith  Eqn . 33.  In  both,  the  variation 
9 E 3 

of  E is  as  l/r"^.  Thus  it  can  be  asserted  that,  in  the 
0 r 

tv;o  dimonsional  situatioii,  the  tv.’o  geometries  considered 
are  just  as  effective  as  any  others  that  may  be  conceived 
for  dlelccl  ropl'iorctic  dc;flection. 


A.  5 . 2 Th roe -Pi meiisi  onal  Crorvtrif-s 

If  a number  of  }')oint  charges  — '^n  located 

at  points  r,  , r^  — r , then  the  net  electric  field  at  each 
^ 1'  2 n 

point  r in  space  due  to  these  point  charges  is  the  Vector 
sum  of  the  field  due  to  each  such  charge. 


E(r) 


Jl 


^ — 
(r 


*n 


- r ) 
n 


From  this. 
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9E 

9r 


47re 


(r  - rj) 


+ . . . + 


■•n 

(r  - r ) 
n 


(r  - r^) 


(A-61) 


Again,  both  E(r)  and  DE/9r  (and  therefore  E 9E/Dr)  are 
significantly  large  only  at  the  singular  points 


r = r 


1' 


.or  r 


n 


Tilso,  since  both 


1 ' n » 
€ -*  0 , 


and 


tend  to  vxiy.o,  the  contribution  due  to  till  other  poi)it 
charges  j.s  insignificant  in  the  neig)iborhood  of  a singularity. 
An  upper  bound  on  E is  the>efoj'e  of  the  forn 


E(r)  - ^2  • (A-G2) 

r ' 

9E/9r  = - 2q/r^  , (A-G3) 


and,  2 

E 9E/9r  = - ^ 


(A-64) 


Comparing  this  result  with  Eqns.  2s  and  35,  it  can  be  con- 
cluded that  there  is  something  to  bo  gained  from  a transition 


froin  two  dimGiioional  to  tlirf-o  diiiicnnional  electrode  confiyu- 
ratione . 


Altlioucjli  the  general  three  diiaensional  geometry  appears 
to  be  preferable,  there  are  a numljer  of  practical  reasonr. 
which  preclude  its  une.  VJhereas  in  the  cylindrical  or  planar 
geometriee  there  is  an  obvious  path  along  v.’hich  particle  must 
be  introduced,  no  sucli  syimnetric  flow  path  exists  in  the  case 
of  a splierical  cliarged  electrode.  Secondly,  because  tlie 
fields  produced  by  tlie  spherical  charge  are  correspondingly 
stronger,  tliey  also  lead  to  brea):down  in  the  interele  trode 
gas  very  easily.  Hence  lower  accelerating  voltages  must  be 
used  to  avoid  breahdo'vn . Thirdly,  j.f  tlie  part.icles  possess 
any  apjrecjable  velocity  v.’hen  they  aie  injected  in  the  vicin- 
ity of  the  electrode,  tliey  \;ill  spend  a very  small  time  in  and 
will  tlnis  be  only  negligibly  accelerated  by  tlic  field.  In 
tlie  t'..’o  dimensi f^na  1 geoirictricj.' , the  paiti  clt.s  can  very  easily 
be.  made,  to  posse  s zero  rad.i.a’’  velocity  by  immersing  them  in 
a gas  .st<.eijj.i  flo  ji.ny  }ia3  .illel  to  the  a.:is. 


l.ic  eauso  of  tlir ; e reason:  , it  i j;  do;;  iiabl to  operate  v-’ith 
two  dimensional  electrode  configurations.  Our  approach  to  the 
design,  o.f  sucli  a system  v;ill  tlius  l'(  confined  to  tlie  two  dimen- 
sional, c y 1 i n u r i c a 1 g eon le t ry . 


A . 7 Design  /ni'proc cli  to  elect  ins _T oparati  on 

The  differential  cquatio.i  of  motion  for  the  cylindrical 
geometry  v;as  found  to  be 

k. 


II  n I 

X + ^'2  ^ ~ ^ i 


(A-65) 


v;here  k. 


Co  (K  - 1) 

,2  4 

P P 


(A-66) 


and  k. 


18n 


pd^  + pdL(a  + 


(A-67) 
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We  must  examine  this  differential  equation  for  the  domain  in 
which  the  particle  motion  is  not  drag  limited,  but  is  instead 
determined  primarily  by  the  electric  field.  That  is,  to  find 
operating  conditions  such  that: 

. k 

X <<  — ±-  . ,, 


Let  3 be  a constant  such  that  3 <<  1,  and  lot  us  substitute 
for  the  previous  inequality  the  equation 


k 


2 


I 

X 


3 


or 


(A-C9) 


(A-70) 


I7e  v;ill  write  this  differential  equation  in  separable  form, 
and  obtain 


3 

k^  X' dx  ^ 3 k^  dt. 
Upon  iiitegration,  v;o  find  that 


(A-71) 


3 kj^  T 


(A-72) 


v;hore  the  particle  motion  occurs  between  and  T.  from  x to 

X 2 1 


^2  * 


If 


(^2  - ^1  ) 


1 1 is  defined  as 


- b> 


Ax 

AT 
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then 


4P  k. 


(A-73) 


is  n -onclilion  for  e5;sontial]y  drag- free  motion  of  the  particle 

V/e  notice  in  the  result  above  that,  since  y.2  < the 

dominant  teria  in  the  determination  of  Ax  is  x^,  the  initial 
position . 

VJo  are  now  able  to  deloi:mine  which  parameters  must  be 
maximized  or  minimized  in  order  to  obtain  effective  dielectric 
separation  wlien  particles  of  different  sizes,  densities,  and 
dielectri  c coirsLant:  a]*c  pres»  iit  in  the  gas  stream. 

From  Equation  73,  we  find  tliat  for  f ield-detcrmiined 
motion  of  the  poii  icJer;,  l.he  ratio  k^/k^  must  ho  as  largo  as 
possihir-.  Thus  this  I'otio 


c o:  - 1)  V ' 

_ n ^ jO_ 

P r/' 


pd^  -i-  pdL(a  i be 


(A-74) 


should  be  :aa::i.ini  zed  . 


before  t'inaiyzing  the  various  piircimet'^rs  indivi.dually , it  is 
preferable  to  cxam.i.no  Uic  cor.-ilrinati  on  of  sorae  of  them  that 
3J-Q  limited  by  othcu  considerations,  lor  instance,  Equation 
74  suggests  that  for  effective  dieiectric  separation,  the 
applied  voltage  should  be  as  high  ns  possible.  However, 
limits  to  the  applicable  voltage  are  imposed  by  the  possibil- 
ity of  eventual  breakdown  in  the  interelectrode  gas.  Thus,  a 
more  global  consideration  must  apply  in  the  determination  of 
the  maximum  usable  voltage. 
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We  will  write  Equation  74  in  slightly  modified  form  to 
enable  breakdown  effects  to  be  analyzed.  From  Equation  32, 
we  see  that  the  field  in  the  interelectrode  space  is  given  by 


E(r) 


and  thus 


(A-32) 


(A-75) 


Since  r is  the  minimum  radius,  the  field  strength  at  the 
o 

inner  electrode  surface  is  a maximum  and  thus  possible  break- 
down will  occur  on  tliis  surface.  VJc  substitute  Equation  75 
into  Equation  74  to  obtain 


E2(r^) 


pd^  4 pdE(a  -I- 
ICn 


(A-76) 


The  various  parameters  that  affect  the  efficiency  of  dielectric 
separation  can  nov;  be  ei^umcrated. 

Since  the  objective  for  effective  separation  is  to 
maximize  the  ratio  Equation  76  v/e  have  the  con- 

ditions 

a)  E(r^)  should  be  as  large  as  possible.  The  limit  to  this 
is  set  by  breakdovm  considerations.  Thus,  in  dry  air  at 
standard  temperature  and  pressure,  the  voltage  gradient 
should  not  exceed  30  KV/cm. 

b)  The  ratio  r /r,  should  be  maximized.  Th^s  suggests  that 

o 1 

the  inner  electrode  should  bo  as  large  as  possible  while 
the' outer  electrode  is  as  small. as  possible. 
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c) 


d) 


ijliould  also  be  raininiizcd  fron  entirely  separate 
considerations.  Tliis  criterion  coincides  with  the 
previous  one. 

The  mean-lree-path  lengtii  of  the  carrier  gas  should  pre- 
ferably bo.  as  large  as  possible.  There  are  two  ways  of 
accomplishing  this  goal.  First,  it  is  possible  to  select 
gases  v;hich  have  large  mean  free  paths.  An  examjile  is 
hydrogen,  with  a mc>an  free  )jath  of  0.118  x 10~^  metres, 
or  1.85  times  tliat  of  air.  Second,  it  is  possible  to 
rarify  tlie  gas  surli  that  a corresponding  increase  in  the 
mean  free  path  length  occurs. 


In  the  analysis  of  die.lectric  separation,  v/c  will  con- 
sider tv.'o  gases,  air  and  hydrogen.  Also,  eacli  of  these  v.’ilD 
bo  consi.dorod  both  at  standard  tempcral.urc  and  prcs;rure,  and 
at  the  reduced  pressvn  e.  of  1/lCO  atmosplierc's . 


Anothr./  desigii  pcs- nc  i i.r , whi  c))  is  Jiot  included  in  the 
previous  discus.sion,  ran  be  evaluated  from  Equation  73.  Since 
the  condi.tion  for  obtaining  drag-free  motion  also  involve.^ 
mi  niniv.ing  the  ratio  /n:/AT  in  tliic  equation,  v.'o  see  that  Ax 
should  be  as  small  as  wc  can  make  it.  . From  tlie  expression  for 
Ax,  v.’r  c'lH  derive  a fift)i  condition  for  effective  dielectric 
separation,  namely 

c)  particle  trajectories  .should  be  confined  as  close  to  the 
vicinity  of  the  inner  electrode  as  possible.  Ideally, 
the  motion  of  tiro  particle  should  terminate  at  the  inner 
electrode . 


These  are  some  of  the  design  criteria  tliat  must  bt.  kept 
in  mind  in  the  actual  design  of  a dielectric  separator. 
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A • C An  Alternative  Approach  to  Dielcctrophoretic  Separati on 


In  an  alternative  approach  to  particle  separation  evolv- 
ing out  of  the  previous  design,  the  dielcctrophoretic,  inward 
foicc  on  the  particle  is  balanced  by  an  outward,  centrifugal 
force,  by  placing  the  particle  in  a circular  orbit.  A suc- 
cessful design  based  on  this  approach  should  eliminate  any 
considerations  of  the  viscosity  of  the  carrier  gas.  We  will 
analyze  tJids  scheme  in  detail  in  this  section. 

We  must  calculate  the  magnitude  of  the  centrifugal  accel- 
eration required  and  consequently  the  angular  speed  required 
to  be  imparted  to  the  particle  in  order  to  achieve  this  dyna- 
mic balance.  An  equally  important  question  is  the  dynamic 
stability  of  this  equilibrium  situation.  Here,  the  magnitude 
of  the  angular  speed  is  derived  for  a given  particle.  The 
problcjn  of  dynamic  stability  is  considered  later. 


A. 8.1  T)ic  Pqui]  ilj'-m'm  Orbit  of  a Paj-ficle 


Thi  c 
circu.l  ar 


ccntripei.al  aiCccO  era  Li  on  of  a particle  moving  in  a 
orbit  of  radius  r v:i  i.h  an  angular  velocity  w (r)  is 


Oj,  - "W  (r)  • r. 


(A-77) 


If  this  particle  is  in  equilibrium  under  the  simultane- 
ous acceleration  due  to  a oielcctroplioretic  force,  this  force 
must  be  equal  in  magnitude  to  the  force  of  centripetral  acceler- 
ation. Tlie  value  of  the  acceleration  due  to  dielectrophoresis 
is  obtained  from  Equation  3G . Here,  the  net  outward  force  F 
will  be  expressed  as  a sum  of  the  (outward)  centrifugal  force 
and  the  (inv/ard)  dielcctrophoretic  force. 


F = m ubr)  . r - m (K  - 1)  ~ 

P P'  r-’ 


In  dynamic  equilibrixim,  that  is,  when  F = 0, 


u (r)  • r = 


(K  - 1) 


.2 


(A-79) 
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or 


This  result  can  be  interpreted  in  two  v’ays. 

a)  The  equilibrimi  radius  of  a given  particle  is  unique, 
for  a given  angular  velocity  profile  u(r),  and  is  deter- 
mined solely  by  the  particle  dielectric  constant  and 
density,  being  independent  of  the  particle  size. 

b)  For  a given  radius  to  which  particles  are  confined,  the 
angular  speed  of  a particle  is  determ-*  by  its  dielec- 
tric constant  and  density. 

Of  these  tv.’o,  the  first  viev/point  is  of  greater  importance 
for  dielectrnphorotic  separatio)i. 

A.  8.2  Sel  f—stabj.lity  or  Stability  in  Vacuirn  of  the  EquiDibrium 
Ofbi~t  of  a Partic3o  ~ “ 

It  is  of  interest  to  d'^ter'"--  vdieiher  a particle  in  a 
dielcetrophoretic  field  in  i nee  of  any  drag  forces  due 
to  a gas  flow  stream  is  in  stable  orbit  or  not. 

In  the  al)sencG  of  a tangential,  field,  there  is  no  torque 
on  the  particle  and  hence  the  angular  momentuin  relation  takes 
the  form 


9 2 

(mr  oj)  = 0 (A-81) 

Assume  a particle  to  enter  the  cylindrical  interelectrode 
region  v;ith  a tangential  component  of  velocity  v , at  the 
radius  r . Then  the  initial  angular  momentum  is  mv  r , and 

^ p p/ 

since  this  must  remain  constant  throughout  the  path  of  the 
particle , 


= mv  r 
P P 


2 , V 

mr  w(r) 


(A-82) 


the  particle,  the  net  radial  force  is  equated  to  the  rate  of 


change  of  radic  ' nomentum.  The  radial  force  has  two  compo- 
nents, one  due  to  the  electric  field,  and  the  other  the  cen- 
trifugal force  on  the  particle. 

Equating  the  sum  of  tliesc  to  the  product  of  mass  and 
radial  accelercition  yields 

-I  1^2 


mr  w"^(r)  - in 


c (K  - 1)  V ‘ 
o o 


as  tlio  equation  of  radial  motion. 


(A-84) 


Several  interesting  conclusions  can  bo  drawn  from  this  result. 


Since 


— - . ^ 

^p 


(K  - 1)  2 

V 

p P'^ 


(A-C5) 


the  particle  trajectory  is  in  general  an  unstable  one,  with 
the  particle  either  spiralling  invaards  v/ith  constantly  in- 
creasing radi.ai  velocity,  or  spiralling  outv;ard  nt  a slov/er 
increasing  radial  rate.  However,  there  is  one  particular  in- 
jection velocity 


V = 

P 


s^(K  - 1) 


pi  J- 


(A-86) 


at  which  the  particle  is  in  a circular  orbit  about  the 
inner  electrode. 

For  an  ensemble  of  particles  which  have  various  values 
for  the  ratio  (K  - 1) , this  offers  a means  of  selection. 
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Particles  with  high  dielectric  constant/density  ratio  will 

while  less  dielectric  particles  will 


spiral  inwards  from  r , 

• , P 

spiral  outwards . 


Rgdi^  of  the  Equilibrium  Orbit  for  a Given  Anaular 
vuiociry  Profile  (j(r)  ^ ° 

It  is  assumed  that  the  particle  is  constrained  to  move 
in  a gas  stream  flowing  in  the  annular  region  between  the  cy- 
lindrical electrodes  (see  Figure  A-4) . In  this  case,  w(r) 
an  externally  determined  parameter  in  the  previous  equation. 

The  orbital  radius  of  a particle  in  equilibrium  is  then  found 
to  be 

Iv? 


- 1) 


V 


oj(r) 


(A-87) 


For  a given  angular  speed  oj(r),  the  equilibrium  radius  r 
is  determined  by  the  particle  dielectric  constant  and  density. 

As  an  exnmp3e,  it  is  assumed  that  w (r)  is  a constant 
(i.o.,  w^)  between  and  This  represents  an  ideal  mo- 

tion of  tlie  gas  streaim  The  dependence  of  r,  the  equilibriu:^ 
radius,  on  K,  the  dielectric  constant,  is 


■’nr 


T 


(A-8P) 


j A.  8. 4 StaMjUty  of  a Particle  in  Equilibrium  Orbit 

distribution  at  equilibrium  will  be  considered  to  be 
|_  determined  by  the  angular  velocity  profile  u)(r).  Each  parti- 

cle comes  to  an  equilibrium  at  a specific  radius,  depending 
j only  on  its  dielectric  constant  and  its  density. 

A given  velocity  profile  u)(r)  can  be  obtained  by  con- 
^ sidering  the  particle  to  be  immersed  in  a fluid  stream  flowing 

in  the  annular  space  between  the  inner  and  outer  electrodes 
(see  Figure  1) . Fluid  streams  of  various  angular  velocity 
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figure  A~<!. 

ANGUIV^R  VELOCITY  PROFILE  OF  7\  GAS  STREAM 
FLOWING  UETVJEEN  CONCENTRIC  FLRCTRODES 
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profiles  can  be  obtained  by  having  a gas  flow  between  sta- 
tionary electrodes,  by  rotating  one  or  both  electrodes,  o- 
by  a combination  of  both  these  methods. 

The  stability  condition  for  the  equilibrium  orbit  of  a 
particle  is  as  follows: 

F(r) 

^^n  first  non-zero  derivative  of  F(r), 

then 

8^  F(r) 

<0  is  a condition  for  stability. 


Here  we  v;in  simply  assume  that 


a F(r) 
8r 


7"  0, 


(A-8r.) 


and  therefore  the  condition  for  stability  is 


3r 


< 0 


Since 


(A- 80) 


J""  2 . 

m~  " 


e V 

o , . . ''o 


----(K-l)  „ 

p P ■ 


(A-Sl) 


3F(r)  { 2, 

— =-  m^ai-(r)  + 2ro){r) 

( 3r 


3 c 


(K  - 1) 


.2 


(The  factor  m is  dropped  from  hero  on  out.) 


(A-92) 
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The  condition  for  stability  is  thus: 


w (r)  + 2rd)(r)  + 

3r 


--  (k  - 1)  —^2“  < 0 • 


(A-93) 

To  make  this  equation  more  tractable,  v;e  substitute  for 
r from  Equation  87,  to  obtain 


w (r)  + 2rw(r)  + 3 (r)  < 0, 


(A-94) 


r 2 u)(r)  < 0 

3r 

on  the  implicit  assumption  that  o(r)  > 0. 

The  procediny  equation  can  be  v;ritten  as  either 


(A-95) 


(i)  _ ^ ‘‘J  ( ^ ) < _ 2(0  (r) 


< - -j- 


3 (0  ( ) . ) 
o)(r) 


v.’hore  o.'(r)  is  constrained  to  > 0, 


(A-9C) 


(A-97) 


oOine  coji clu.s ions  tliat  can  dra'.vn  from  i-hese  cquatioju’ 

are : 

a)  A uniform  angular  velocity  profile  (vs.  radial  distance) 
is  not  sufficient  to  obtain  dynamic  stability.  Since 
3co/oj  is  then  zero.  Thus  the  idealized  case  of  uniform 
angular  flow  considered  previously  is  not  sufficient  for 
orbital  stability. 

b)  To  obtain  stable  particle  orbits,  3(o/(o  must  bo  nogati.vo. 
Since  co(r)  is  positive  by  implicit  assumption,  it  follows 
that  3(i)/3r  must  be  negative.  Thus  a configuration  which 
produces  high  angular  speeds  toward  the  center  is  required 
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if  stability  is  to  be  achieved. 


Such  a configuration  could  be  obtained,  for  instance,  if 
the  inner  electrode  v/ere  rotated  with  respect  to  the  outer 
electrode . 

A.  8. 5 Possible  .stability  Conditions  for  a Specific  F'low  Profile 

It  is  possible  to  solve  the  preceding  ecuation  for  any 
specific  radial  flov;  profile  u(r)  and  thus  to  predict  tlie 
rcgiojis  in  v/hich  a particle,  in  a circular  orbit,  can  be  ex- 
pected to  bo  st^lble . As  an  example,  assume  thc'it  a flov;  profile 
is  created  in  which  the  gas  stream  flow  velocity  as  a function 
of  radial  position  is  expressed  as 


w (r) 


(A-98) 


This  velocity  profile  is  obtainable,  for  instance,  if 
the  outer  electrode  is  stationary  v.-hilc  the  inner  electrode 
spins  wit.h  an  avigular  speed  of  giving  rise  to  a linearly 
decrcasijig  velocity  profile  as  a function  of  r. 

Using  Ko\iation  96,  the  stability  co.ndition  is 


cir  o 


r.c. 


(I'l  - 

_2_ 

ii\ 

r) 

<'-l  - 

r 

lU 

o 

^^1  - 

r ! 
o 

1 - < - 2 

“o 

- r) 

i - 

- 

(A-99) 


(A-lOO) 


This  simplifies  to  the  condition 


p (r^  - r)  < 1 


(A-101) 
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J 1 * (A-10?) 

for  ntaolc  orbiV.s.  That  i.'j,  tlie  region  of  ntability  for  a 

linear  gan  flov;  profile  ir  tlie  outer  one-thircl  region  betv;oen 
the  tv.’o  conductors. 

The  above  dcriveo  .stability  condit.ion  is  .son'.cv.'li-at  sur- 
prising in  one  respect;  namely  that  the  .stability  of  the 

is  indepcr.c;cn t of  ^qI  the  miagnitudc  of  tlie  annular 
speed  v.’itl)  v.’hich  the  inner  electrode  is  rotated.  In  theory, 
stable  particle  motion  could  be  obtained  witli  very  lov;  rota- 
tion speeds  of  this  inner  electrode.  In  practice  Jrov.’cver,  a 
ininiiuum  angular  rpoco  must  Ijc  vi5;ed,  Tiris  is  true  be- 
cause the  .stabil.ity  condition  v;e  have  derived  is  one  of  the 
bounc.c(;-- input  Irounced-output  kind,  and  says  nothing  cil:)OV^t  tJie 
cxcursion.s  of  the  particle  about  i t.s  .stable  o.rbit.  As  the 
angular  speed  is  dec.-r'na.Sv.d , tlie  oxcur.sions  of  the  prarticlo 
around  its  equilibrium  orbit  become  larger  and  larger,  and  at 
some;  poi.nc  the  uncertainty  in  the  particle  orlrit  booc)i,ics  com- 
*->aranle  to  the  tot<.il  v/icitli  or  the  stability  region.  mini- 
mum angular  .speed  for  the  inner  cloctrodo.  must  therefore  ' 
be  used.  Til  is  minimum  veilue  of  iau.st  be  di  rived  from  sepa- 
rate, perhaps  csyporimental  consideratio.ns . 
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A . 8 . 6 Improved  Particle  Stability  Conditions  for  the  Flow 
Obtnined  by  a Rotation  of  the  K.lectrodes 

It  is  possible  to  solve  the  stability  condition  derived 
in  the  previous  section  for  any  specific  flow  profile  oj(r), 
and  thus  predict  the  existence  of  regions  of  gas  flow  in  which 
a circularly  orbiting  particle  can  be  expected  to  be  stable. 

In  this  section  v;e  will  consider  the  gas  flow  profiles 
resulting  from  rotating  the  cylindrical  electrodes  about  their 
axes,  and  der.ive  the  corresponding  stability  regions  for  the 
particles.  The  gas  in  the  annular  region  is  assumed  to  have 
no  externally  impressed  velocity  profile.  Thus  the  angular 
velocity  of  the  interelectrode  gas  is  determined  solely  by  the 
speeds  of  the  inner  and  outer  electrodes,  which  are  assumed  to 
be  c'-ncentric. 

Such  a gas  f]ov;  must  necessarily  be  3 anv’ nar  and  satirn'y 
the  governing  fluid  dynamic  equations  and  boundary  conditions. 
However,  these  conditions  are  not  sufficient  to  ensure  tlie 
stability  of  t.he  gas  flow  itself.  The  instability  of  tlie 
flov;  of  a viscous  gas  bet'/cen  concentric  cylinders  is  a kind 
of  inortxal  instability,  and  its  possibility  must  be  taken  into 
account.  This  will  be  dealt  with  later  in  the  analysis. 

In  tlie  analysis,  we  v;ill  consider  that  the  inner  and  outer 
electrodes  spin  about  their  common  axes  with  angular  velocities 
Wq  and  respectively.  Under  sucli  conditions,  and  assuming 
laminar  behavior,  the  velocity  profile  of  the  gas  in  the  annular 
region  can  be  shov/n*  to  have  the  time-average  form 

* Taylor,  G.I.,  1923,  Phii.  Trans.  Roy.  Soc.,  London,  A223, 
p.  289-343  
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w(r)  = A + 


(A-103) 


where 


2 ■ 2 
"^0 


(A-104) 


2 2 
“^0 


(A-105) 


We  v;ill  npp.ly  tlio  previously  derived  particle  st^-ljility  con- 
dition, namely 


9 cj ( r ) ^ 

8r  ' r 


(A-106) 


to  this  particular  veDocily  profii.e  and  obtain  the  re.cjions  oi 
stability  ioi  partic''e  i..otion. 

f'lOD.viiv^  io.v  boUi  ; ;idcr;  of  this  cc’\'Oti<  n i.>nd  (iqiu  lincj , 


we  have 


2J^  ^ n h < \LI 
- 73  - - • ■ ,3_ 


(A-107) 


as  the  stability  ccmd'tion;  t)iat.  is. 


- < 0 
r 


(A-lOO) 


Since,  in  the  expression  for  A,  r,  is  always  greater  than 
2 

Tq  , the  condition  for  stabilrty  thus  reduces  to  the  very 
simple  one; 


(A-109) 


, “0*^0 
“1  2 


(A-110) 
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Wc  will  for  siinplicity  suppose  that  the  outer  electrode 
is  held  stationary  = 0)  while  the  inner  electrode  is  spun 

at  the  speed  Wq.  This  situation  satisfies  the  stability 
criterion,  and  the  particle  orbits  are  then  stable  at  any 
radial  position. 

A . 8 . 6 Inertial  Stability  Conditions  f or  the  Gas  F low  Profile 

In  order  to  obtain  slab.le  orbital  motion  of  the  particles, 
it  is  necessary  to  rotate  the  inner  electrode  at  an  angular 
speed  (Oq.  Although  tliis  arrangement  stabilises  the  particle 
motion,  v;e  must  additionally  ensure  tliat  tlit:  resulting  gas 
motion  is  itself  stable. 

As  dir.cussed  previously,  it  is  possible  for  a gas  stveam 
under  such  boundary  conditions  to  exhibit  a form  of  inertial 
instability  arising  from  the  accelfe:ation  of  the  container. 

The  follov;ing  is  a condition sufficient  to  ensure,  the  abseric;e 
of  thif.  kind  of  si. ability J;nc-v.’a  rr  Tiiylor'  instability. 

For  stable  gas  flov;,  v.’c  must  have 

~ r^co(r)  > 0 (A-111) 

dr 

in  the  entire  region  r^  < r •$  r^. 

V7c  substitute  the  knov;n  form  of  w(r)  from  tlie  previous 
section  to  obtain  the  condition 

(Ar^  + B)  > 0 (A-112) 

dr 

which  reduces  to  the  condition 

^ ^ ° • (A-113) 


* C-S.  Yih,  "Fluid  Mechanics",  Ch.  9,  Sec.  3,  McGraw-Hill,  1969 


Block 
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I _ This  condition  for  the  stability  of  tho  gas  flow  is  the 

; I exact  opposite  of  the  condition  required  for  the  stability  of 

) the  particle  orbits,  derived  in  the  previous  section.  V7e  are 

: ^ thus  forced  to  conclude  that  the  conditions  for  stable  gas 

j flow  and  stable  particle  orbits  are  mutually  incompatible. 

I 


J I * 

I 

I 

I 

I 


Reproduced  from 
best  available  copy. 


Apivcndix  3.  FQRTR7iN  Pr 
in  the  Pla 


rograms  to  Genor^^to  Particlo  Tra joctoric 
^nar  and  Cylindrical  GcoinotriG.s  ~ 


Tf5":  riixL' 

nr/i  KAr.”A 


rC’C.SSK-  12 
r.:o-  :i3 
KAP.''A=IC.e 

PP=  < )..0Z(  ?.r*50.C-l.J))»«2 
ZT/.-C,  1S?»A 


1 

poe:j=  C.9 

ACCEPT  ' .VTAT  CTCPJ;  3.0  FOn  IirFAVET)  ",  ETA  1 

If  <;:ta!.!  1.C.1;)  et,-  = etai 

f 

! 

ACCENT  "lor*;?  ",  rcr:.'i 

IF  < r 0 c.  1 . f. n FOf'.'=iof*ii 

} 

C 

WF.I  1 F(  ID,  13) 

10 

FOl.'lATC'  !;.LCir.0CE  5P/.C1NG  IH  Cl.  7 CFA.S) 

) eauc  1 1,  roil 

S3 

rOm.'TCF.l.J) 
M»H/ ! IT. 0 

vri  t;:<  1 r,  05} 

05 

roivATc  r.'.pfiar.  D:r;f.Tr..'i  ir:  i:icr.:;;s  7 crs. 
rr/c:  11,  30  53 

33 

roiM,-a..r5. 1) 
i't=  < rr  • 1 5-6) « } 0 
U!  I :r.<  1.',  35; 

35 

rCllC/TC'  VCLT/'F  7 CFV.DV) 

. n-,'.;  c 1 1;  .'..')\o;.'is 

n « 

63 

FO.'.ATC  • 7.  1) 

C 

63. ' K/'  ? . f • t ( )',■.'■  r,'  - ! > •'”31.7  f-  0 6!'.  10/ VS 

C 

u.'.T/,=  I'.; . c« : f/rriio/cs 

• 1 3.':.'  = ro.';:'i;.  1,  c.c 
va,=  c,3  ' 

ACC::= : -ir  • ' : £;;/  c 1 . c-:  3 £r;*-i  0 •*  3 
T/CC'=;  cc:.‘H=  IC3.  e 

c 

v.'.iTrc'ic,  VO 

• 

53 

1 

roJi'./.T ( : 3,  • ; . r I r; ci;  ’,  TIC,  ”.’.'i.rci tc?,  'ace: 
Tat,  •>  1 ?;.L’  Aca. ',  vi  5, ’riMAC' acc;;  tiiiev  ' 

0 

T3,  • c:-!.  ',  7 IC , 'C'./  5U  ST:  ),  •CI/5EC50',  ' 

3 

T'.r,  • C'.VCECVT’, i vv,  'E::/5;.c5o  ;:.5ecv> 

vr : Tr.c  I.',  f c) :: : r.;,  T'.T!.,  7,7Cf)  ' . 

6C 

} C'MCaTc  7 :i,  E 1 1 . 6,  7 1 5,  L 1 1 . 4,  T 0;',  E 1 1 . 4,  TA  1,  1 1 . 1 

65 

70 

C 


Tj 

!;r=  c.  0 

LO  Vi'  1 = 1,  ICB 
Ti;ii>7i:;[*p.  i 
ro  65  j=  I,  ice 

ro  5U  ■ I 1 1.-  6 • VET.,*  0,  5*  < 1 E~6)»»  ACOJ 
Vl.I.  = Vil.«<  Ii:-6)»AC&'1 
FACc;;  = - AUl  ,'(/yi  a 

DAcr;j--r5TA»VLi. 

ACcii=rACc:j-  i.Atcv 
Tro£N=;'0.‘:iMi*  it'3.0  ■ 

•<■^•(.1.=  v;'L*i(»  lii’.o 

I AC(,M  = ACc.'J  = l!»  IP.T.O 

T>'Acrj  = F,',cc:;»ii.  icc.  o 
“CAT.';-  r.ACcmi*  iDB.  o 
Cc;jvi:;cE 

VF.I  TEC  ID,  63}  TPOS,',  TVU,/  T6,CCN/  IFACCit  TDACCW,  TIME 
COM T I HUE 


STOP 

DIO 


I 


TTPE  rm.D3 

REri-  KAPPA 

c . . 

EO-B.aSE-  12 

KHO=  3t3 

KAPPA-' 10.  0 

PP>  ( AEOG(l>'J.fl>)»»S 

ETA’C.  Ill-iE-  A 

por.-i^c.  1 . . I 

ACCEPT  "VTA?  {TT-PF  C.O  FOR  t'EFAU.T)  *’»  ETAI  ] 

IF  <;.TAl.t;E..''.C)  ETA-F.TAl 
ACCEPT  "F0^:.•7  "/P05U1 
IF  {POSN 1 e. 0)  roc.-i'-ror?) I 

c 

vnnrt  m,  lo 

10  FOIMAK'  EXECTr.ODE  ri/.-lETEE  111  Cl.  1 (F4.  2)  •/> 

11EAD(  1 1<  Of  )H 
£0  FOniAK  FA.  2) 

1I  = H/P/J  Ol’.  0 
Vr.IThC  10.  £5) 

25  FOIMAT.-'  -PAPTICEE  DI/.-lETER  Itl  111  CHOUS  7 (PS.nV) 

READ(  1 i;  30)  C2 
30  FOlMATCFl.l) 

CO'-  (C.'t  IE-A)**P 

VRITECIO.  OS)  . I 

35  FOI.'IAIC  VOLTACF  7 (F7.1)V> 

Rl.AtX  1 1.  AO)  VOLTS 
AO  FOI:iAT<K7.  1 ) 

c 

ALPHA'Ei'  ClAPPA-  l)*V01.TS»*2/f.l0/P2 
ALPHA-  /!.PM,VH**  A 
LET  A- lU.  V' f.T,Vr.Ry/H2 

c 

Tro:N  = por!!»K*  ico. o 

TVF.L“C.C 

VIX’O.O 

ACCN>’-r,i'ii,vros-.*i3 
TACCP-'ACC  ;•;<»  ica.o 
c 

VPITEUC.  SO) 

50  F0F:1AT<T3.  '!  OS!  TIC!)'.  TIA,  • VIEOCI  TO' , T.'-'J,  ’ ACCEiEIlATl  ON”. 

1 TAP.  ’FIH.n  ACCN-'.TSE,  'PPAO'/CC;!  TJMFV  ' 

2 TO,  • 0-1.  •,  TIC,  ’ C;/rEC.  TO).  • C.1/5ECS0*. 

3 TA2.  •C.'!/;-ECSO'.  TSL,  'fcVSl  CSC  MSECV) 

vf.iteuo,  co)'iPoy:.TVi:..7Acci 

CO  For:-!  AT  < TP.  1,1  1 .A.  T IS.  El  1.  A.  7L3,  El  1.  A.  TA 1 . E 1 1 . A.  T5  A.  E 1 1 . A.  F7.  3) 

C 


T, 

*nE«c.o 
DO  7C  1= 1 . 1 CO 
T1I1L''T1!1F«0.  1 
DO  65  1.  ICO 

P0  5TJ=1’0':;*  lC-6«VF,L*0.5tt  1 E-6)»*2»ACC71 
VlL=VIT.*<  ll.-C)»ACai 

FACOH-/J,rT!iVrOS.’J«*3  I ■ 

DACCt.'-'C-F  1 AiVLL  I 

ACCtl-'FACCi- PACCN  S 

Tl’OSN''-POSM.tl«  ICO.O  ’ 

TVFl."  VELAll.  1 CO.  0 I 

TACC'J"  ACt'.ri»ll«  100.0 
TFACC'.-  = FACC!1*II*  100.0 
TDACC.H-'DACC:;»ll»  100.0 
65  COHTINLE 

WRITFC  10.  6!1)TP0EN.TVEL.  TaCCN.  TFACW.  TDACCN.  TIME 

70  C0.NT1MUE  ) 

' STOP  i 

END  t 
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I 

I 

I 

1 

I 

I 

i 

I 

i; 

f: 

i: 

i: 


i. 


! 


APPENDIX  C , FORTRAN  Programs  to  Generate  Particle  Trajectories 
and  betermine  Particle  Distributions  as  a Function 
5^  Particle  Characteristics. 


r.Tt 


c 


#ri 


«#!:r 


Affrr?  <|)  er  rLfr^nry  rc  •» 

j.rr.r;  'vv7:r  ::.»c7rorr  r/.rirr  n:.  ? 

t “ir.r.T!  :!.';rv:\jr!  itAfirr  r;  r-.  •* 

!■<«•  virr  ■ 

ik:;,  .r.-.  i>r*.-_cr<-i/rr) 

IKV.r  tr.r  I, ;■«.  ,,fj 

Arc:.-?  iojjtic;:  ik  c:.  ? 

. Of:’*Tr^«:;/lir/:  « 

ACfiTT  kT  !;rTrjcs:<s>  ? 

I!  r:.  i;  r;.-.. 

/crrt’7  •*rrj!?5v“7  p*  prioFniyr.r'  ? 

I fc :*r# rr*  i )?n.*: r*  ri  't;* 

! j r r>r ~ . 1 1 . r- ( ».  r.: 

•*^c.-y7  •TATTirr.r  ci.';iE7rp,  iri  ^ 

AfriTT  "rir.rrtric  rc:;r-.;:T  ? 

sc-i  tr.;si-:r  ix  c.:/c:-»?  ? 

Acc-rt  “AF?:.i2r  voltags  ? 

E2=f.5t--is 

/ip:/.«rc«<A."Ar?;.- 1.  r>«v«es/r5.:o/r.j« 

!rtr..;.>7/r::..r7.  j£g  •tg  c 

r.:‘  c-  a/.:.  : > 

* .*T«*  . • » 

TVi’.-t.r 

irc;:;..  r7.  s)”i<  c?0KJ-r. 

/7f  . = 

7^.CC. ACC £. !*'t  Is* 

vr.ir?.{  ir.  ii> 

v3r::A.T(/////7;,  •:-3.'iT!c:  •.Tit,  •v'*  r-t -.r  

TAi,p:::.r  rcc^’.r^J’  ' 

-T'.’V;*  •'■7:-  t /jrcr:-, 

'•.!7E<i:.ic!7-;r  ,r  T/xc:'  ---t  /j 

ys.cA.-t-i,!!;.  *■  I-,:  11.  A,  7r3.ru. -I  n.  A, 

;:7i;;r»: 

:i  c:  jri.itr 

to  Z'O  #?«  !#  j rr 

ir-c  ''.n.+ s"-  iF.-rcr:: 

• rj.*"  2-*!z-c^/cc: 

irr::.,z:*  !.  :/:c 

*-'**>- ? 

t-ir*-::?/."’ 

Att; 

I 

77/:  Dca  7o  ac 
7:-7;  >r;r:>ri<- 1>? 

v.f3,»v!i.»5..|,,|rr 


c 

Cl 


l.-tiy.-.cc  *71  irr 
7:'Af (..’■•r acc:.'i»7i  irr 
• rAtc:»t'cc’  »"i-  ir? 

7!:.i":;7;;:-  .|:--2 

•s*rs»7r  o;:’.TV*i.,  Tacoi.tfacc:;.  Tr./.coi.  ri-'z 
ito'jru.LS.rt/ritcc  7o  te 
a;.Ti::ix 

:tor 

ffc 


C-1 


I 

i. 


!T  ri!Lrr 


r.rriTT  “tru::!  r tf  i>  or  n.,-:;/.;  tr>  rr;;;:T!*-  ? 
,-cc-rr  “ot~T  :!.>i'"rrr  r.-rif:  ir  ? 
ri-n/irr 

/cc:rT  “If.;-;  •LUTi.'ri  r/ji;r  r . ? 
rr*:  .virr 

l>:'>-.".nr{ri/rr) 

;i/rr-  i.rs 

fccr7  roMTi.->:'  ir  r:.  : 

:e!:.-''.ror:/i:r/ri 
/'CcirT  •v.irci)  or  •:'‘vrorr:r<r>  ? 
ire:!."''.  !'rr-.*. 

IFC;:: .rr.r>r'.'-  'rTi-.; 

ACriTT  ir  /~;cr::!:r.rr  ? 

s?rr*i,*‘  '•  r** 

iFc::'-.  rr. 

.'.rcvTT  "r/r.Tif;.:.  i;:  rirr/"  ' ■» 

r«t'!'-c 

A.rf:  ?T  "ri  ?i.ycT:-i  e rrr  ?■.•;.■?  ? 
rrcrr-T  "pr.-nr-r  rrrfiTY  i:;  r:/c:-'»  ? 
f;-;c«?::o*  I • r 

A.ccrrT  voriYCE  ? 

iC»f..r:£-i£ 

/LrH,-«£r»</;:f.rr;>-  i.  r>*vi»r/?;:a/rv--  a 
XT.r 

>Ftr.A.‘^r/:-:i..rT.  >t'.:)co  tc  i 

F,f  If  .rn;r//cn:C‘-  r‘■«£»r:^cn^■-r3^  -tr.AfTr.Sfv::)) 

vn.=c.r 

TVlir-C.? 

IFcrA.i'.  ij;> !.  r/p:t::«’r3 

I r«:’A. !. r- r? r;-) /< "o for: 

VACC  'ACOnir  !•  f 


vriTrf  ir,  i:> 

rcnrATCA/T:.  •?cn?!3;".':ir.  •••r!,3c:r 

TAi.’Fi:.!,r  ACC::*#  •;£/;.  • ;-A3  .-rr;  T!--v 

T3,  ’C:.  •*'.  If,  'rvf re 'c 
TAi,*cvr:cf rvrrff:  ;;Jrcv> 

i.Ti  Tr<  I i,  £c>  7: 3 rir!..  7/cc; 

KOflrf  TC7£,£I  I.  A,  711,  El  I . A,  T?i,  Z 1 1 . 7«  I,  £ I I.  A,  7£A,  £I  I.A,  77.  J> 

i:*l 

KTirrtf 
GO  7C  A£ 

:i*r*  n 

vofr;«ro£r!4  ir- «*••:-£♦£.  ez-  i s^-acc; 

Vfl.«Va.»l£-tA  ACC” 
lrc:AA.£7.  !>:>  i.r/!'7f:>»8 

!.7-rct!.)/iror:.-c.£i:?7rt;».£>j»j 

rrcc:*-fz?;-:A.,:: 

CACC-”-ZZTA*VEI, 

Acc‘-r/.cc:»r/,cc:: 

IMr0ri:.CT.'£/7I.A;;£.KTIKl.lI£.-:>73  73  «J 
7-0£;>Por’:»M»iz:; 

7viL»va.»ri»  iFf 
7/.CC>ACC;»7l»  170 
TE/.corA.cc  vr.i,  in 
7iiArc!»rACC>ri-=  IE? 

7Ii:e-kti:ie»iz-:- 

•,Tirr.<is,rc>*.!c::i,T'  ZL,TAecii,irACc:-,TrAce::,Ti:;: 
lt<ror!;.cT.;c/M>co  lo  ae 


C-2 


ri.-.u>« I,'.?  . 

rjA«f:*f.rr 
riAc  ic 

/rrrrr  or  n.mr~tr>  r;c::iTrv  » 

«r:-'.~T  "oiTir.  r;.rcTrcsr.  ’•/.riur  i;.  c:.  i 
ri*n/K: 

/rt!-r  “::;:;rs  a.:- errors  taiiv?  i:;  c:.  t 
r.r«r.:/irr 

ir<:-:.'.rr.  i>r»«.orcri/;  :> 
irc;:/.ir.  r>r»rt.?c(r*r  i/rc-  i.o) 

/.fcrrr  "imti/x  rrriT:;::  i;;  c:.  ? 

7ro!:.-«At-o':: 

T05:>rror;/isE/r,i 

Acerrr  ■■,•!'<!>  or.  Kvrrrcsr'csi  i 

ir<::r.rc.  i>r7A«r.  i;-rr-A 

ir«:::.r:.E>STA-f  .t 

Arc;??  “rrr.rji'r.;;  jr  Ar::orrHEr.E5  ? 

ircsF.so.  i>?n.«:.r;.*r--;si,*.r;: 

ir<::r.rr.r>rn.«r.  i !5"-{-»a~: 

Acerrr  •TAsrie.r  r'-rsirv  i;:  ci/cie--:  ? 
r;-;:«?::o*ir3 

ACCirr  ~A??LIEC  VCS7ACi  » 

Er*s.cjt- 12 

EC  C I>i.2 

A;:Arr/>!r!E!-ti> 

7ro5:;eAT"cr:: 

P05::*=r?cr:;/iE2/n 

Ar.?::A«EXA”APrA-  i.r)«v-*2/ir:o/r f/riu.^4 
S'l.: 

ir<r.«A«r/rr..rT.  ij:.c,)co  ro  s 
j:'n;?c-:.44i'r/rrL> 

ri*.'A» >•  ••■■'£tr-r'o»E*FFr«{r./.f,*r.3i»x> > 

v;-x*:.r 

r.FLi^r.r 
ircKA.FC.  !J”e 

1 1.  re  :frrr;.=  »rj'>2 

TACC>/Cr;>F.!»  IF£ 
vrirFtir. let 

ro:::/rf//r:.‘rer!rjr::',rii,  "•a.cc'rv.r:;,  •,-cr.  •. 
TAi, ‘rK-xt /re:' •»  rr/„  • t'.-c  -tr:  vtv 

T?, *c:.  •.ri:.  rrr,  •:  /rrer' 

*/.!. ‘rs/rrer' • c;/:;c:'  ".efcvf 
i;,fr>Tr-c:",-.-  ;r..v.Ts: 

ro'v:/r  t iE»  E 1 1 . • I r 1 1 . A,  rre*  E I! . 7/ :» r 1 1. «,  r:  a, 
te  5:  zvi:;E»i.«iir 


P«!!>r7S:.  l>.-C»Vri.*-.  .15-  12-ACCI 
v:r.»vi-.«  ir.-t./ct;: 
ijcr/.r'.  !>•;»  i.r/rer  rt 

I! CEA.i'. 

FACc;«-/'r.:::A«:: 

r.A7f.>-:E7A»vr:, 

ACcr-r/Arc/cc" 

Ti:-;.:;ri:::-ir.-i 
7:c:>:-cr;>M»  i»{ 
v.i!.*-":".*’  i-irr 

7/ti;  •.•f.r.:  ■-i.'.r: 

ir/.ct  •‘F/te;t»M<  ir? 

Tl/tr 

viii:  ii:  ,r'»7rAs:  ruiAtci.e.-fcr :.7:Aec:,7i::F 

MiircTiie  lA.jT.  • •.FA.r-, • PArtia.r.  riA  • vrs.»//i 

mr.r 

eo:i7i::< 


